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ABSTRACT: Treating the wind-excited motion multistory building random 
vibration problem, the input-output spectral density relationships are derived using 
transfer matrix formulation. this analysis the motion the building restricted 
the direction parallel the average wind flow, the along-wind direction; and the 
response represented displacement variable and force variable 
each story unit. 40-story building used numerical example. 
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ABSTRACT: The change the internal pressure building caused the sudden 
opening window door under the action wind analyzed using the Bernoulli 
equation fluid mechanics. The problem treated several levels sophistication, 
from the simplest quasi-steady isothermal incompressible solution the most complex 
solution assuming unsteady isentropic compressible flow with due regards the inertia 
effect. Following the break window door wind, the internal pressure reaches 
the external pressure the opening rather rapidly. After the internal pressure has 
reached the external pressure, there overshoot and subsequent oscillation the 
internal pressure, spite the constancy the external pressure under steady 
wind. The results are compared solutions obtained from the Helmholtz resonator 
model acoustics. 


REFERENCE: Liu, Henry, and Saathoff, Patrick J., “Building Internal Pressure: 
Sudden Change,” Journal the Engineering Mechanics Division, ASCE, Vol. 107, 
No. EM2, Proc. Paper 16172, April, 1981, pp. 309-321 


16192 TORTIONAL VIBRATION EXCITED STRUCTURES 


KEY WORDS: Aerodynamic forces; Aeroelasticity; Dynamic response; 
Random vibration; Random vibration analysis; Structural dynamics; Torsion; 
Torsional vibration; Wind direction; Wind forces 


ABSTRACT: method for analyzing the torsional vibration along-wind excited 
structures presented. The method based random vibration concepts and yields 
the expected maximum translational and torsional responses single-mass structure. 
The approach similar the gust factor method. Aerodynamic admittance functions 
are derived which are used estimate the spectral density function the random 
torque and the cross-spectral density function the force and torque acting the 
structure. These are required for estimating the structure’s translational and torsional 
mean square responses. Results for several examples indicate that the dynamic 
torsional response increases follows: the width the structure’s exposed face 
increases, the structural geometric eccentricity increases, and torsional natural 
frequency decreases. 


REFERENCE: Foutch, A., and Safak, E., “Tortional Vibration Along-Wind 
Excited Structures,” Journal the Engineering Mechanics Division, ASCE, Vol. 107, 
No. EM2, Proc. Paper 16192, April, 1981, pp. 323-337 
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KEY WORDS: Bridges (girder); Curved beams; Finite element method; 
Inertia; Natural frequency; Stiffness; Structural 
engineering; Torsion; Vibration 


ABSTRACT: general finite element displacement formulation presented 
determine the natural frequencies spatial thin-walled curved girders including the 
warping contribution. Also included are the rotatory inertia effects with respect 
flexure and torsion, and the effects antisymmetry cross section (found 
significant). The stiffness and inertia properties the finite element are obtained using 
the solutions homogeneous differential equations governing the static problem 
deformation modes. The obtained element relationships were programmed for use 
digital computers, and few example problems were analyzed. Some charts are given 
expedite the determination natural frequencies horizontally curved girders. 


REFERENCE: Yoo, Chai Hong, and Fehrenbach, Jon P., “Natural Frequencies 
Curved Girders,” Journal the Engineering Mechanics Division, ASCE, Vol. 107, 
No. EM2, Proc. Paper 16177, April, 1981, pp. 339-354 


16193 LOGICAL ANALYSIS STRUCTURAL FAILURE 


KEY WORDS: Bridge failures; Failure (mechanics); Human factors; Logical 
elements; Logic design; Random error; Safety; Safety design; Structural 
behavior 


ABSTRACT: engineering only the product, the hardware, physical system; the 
system which designs it, produces and uses human and therefore complex and 
vulnerable. Human error one form another major cause failures. One way 
attempting reduce this find ways distilling conclusions from past 
performance logical analysis available information. This analysis may then 
stored computer data structure and used basis for predictions regarding future 
projects. The vehicle for this process, the handling rather imprecise information, 
fuzzy logic. logical hierarchy simple deductions has been set represent the 
three main influences structural safety—human error, random hazards well the 
system uncertainty reliability theory. The hierarchy used obtain measures 
safety for two famous bridge failures. 


REFERENCE: Blockley, David I., “Logical Analysis Structural Failure,” Journal 
the Engineering Mechanics Division, ASCE, Vol. 107, No. EM2, Proc. Paper 16193, 
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16200 ULTIMATE LOAD BEHAVIOR CURVED I-BEAMS 


KEY WORDS: Curved beams; Deflection; beams; Load bearing capacity; 
Loading tests; Load limits; Mathematical models; Matrices (mathematics); 
Plastics; Stability; Tests 


ABSTRACT: method analyzing the large-deflection behavior the horizontally 
curved I-beam, both elastic and inelastic regions, using the transfer matrix 
method presented. Three kinds separate loading (e.g., equal end moments, 
concentrated load the midspan, and uniformly distributed load) are considered for 
this analysis. The ultimate strength determined for several given curvatures. 
Experimental studies are also made for six supported I-beams determine the effect 
curvature the behavior ultimate load-deformation. Finally, formula also 
deduced predict the approximate ultimate strength the curved I-beams. 


REFERENCE: Fukumoto, Yuhshi, and Nishida, Susumu, “Ultimate Load Behavior 
Curved Journal the Engineering Mechanics Division, ASCE, Vol. 107, 
No. EM2, Proc. Paper 16200, April, 1981, pp. 367-385 
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KEY WORDS: Elasticity (mechanical); Fatigue (materials); Fracture 
mechanics; Fractures (materials); Fracture strength; Fretting; Mathematical 
models; Stress (mechanics) 


ABSTRACT: The distribution stresses typical two-dimensional fretting-fatigue 
problem has been analyzed assuming homogeneous, isotropic and linearly elastic 
material. “Closed type” solutions have been obtained that are “exact” within the 
assumptions used. According the numerical results presented, small region the 
plate immediately outside the contact area would appear most critical from the 
viewpoint crack initiation. However, under certain conditions (e.g., the presence 
micropores other types discontinuities) cracks could also initiate some 
subsurface sites. 


REFERENCE: Chung, Jung S., Oran, Cenap, and Hoeppner, David W., “Elastic 
Stresses Fretting Fatigue,” Journal the Engineering Mechanics Division, ASCE, 
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KEY WORDS: Anisotropy; Composite structures; Dynamic structural 
analysis; Fourier series; Laminates; Orthotropism; Plates (structural 
members); Sandwich laminates; Stiffness; Vibration 


ABSTRACT: The modified stiffness method unbalanced anisotropic sandwich plates 
extended and applied the dynamic analysis sandwich plates consisting 
core and two equal thickness anisotropic face plates. The formulation 
allows for the dynamic analysis certain types sandwich plates that were not 
amenable rigorous solutions due the existence anisotropic coupling stiffnesses, 
such those with inbalanced angle-ply faces. The method independent the 
boundary conditions and can combined with variety analysis techniques. Using 
series solution, results are generated for simply-supported sandwich plates with 
unbalanced cross-ply and angle-ply face plates. addition, the effect coupling 
the static and dynamic responses presented. Coupling can influence the structural 
response significantly. 


REFERENCE: Ibrahim, Ibrahim M., Farah, Anis, and Rizk, Nabil F., “Dynamic 
Analysis Unbalanced Anisotropic Sandwich Plates,” Journal the Engineering 
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ALONG-WIND MULTISTORY BUILDING 


INTRODUCTION 


Wind excited vibrations building can cause failures structural members, 
discomfort its occupants, unsightly cracks the interior claddings, etc.; 
therefore, important provide adequate design minimize the chance 
excessive stresses and excessive displacements. Since all natural strong winds 
are turbulent with randomly fluctuating wind velocities, probabilistic approach 
most suitable for the analysis such problems. The pioneer work Davenport 
(4) made use the general methodology random vibration which was becoming 
increasingly more popular during the 1950s and 1960s primarily aeronautical 
engineering applications. Since then, numerous papers have been published 
the subject (e.g., Refs. 9-12). Typically these previous works, the motion 
building expressed terms normal modes, the obvious advantage 
being that the solution complicated multi-mode system immediate 
extension that single-degree-of-freedom system. However, several practical 
shortcomings remain, including the following. 

For complicated structure, the determination normal modes can difficult 
extremely time consuming; therefore, the true advantage normal mode 
formulation depends whether not adequate represent the structural 
response only the first very few modes. This implies that the normal mode 
approach may suitable for the analysis wind-induced displacements, but 
may not for wind-induced internal forces. The forces are obtained 
from differentiation (or its equivalent) the displacements, thus higher accuracy 
required for their computation. Secondly, the normal modes are, definition, 
associated with undamped system. The magnitude damping each mode 


must assumed rather arbitrarily lumping together various contributions 

Dept. Civ., Mech. and Environmental Engrg., George Washington Univ., 
Washington, D.C. 20052. 

Dept. Aeronautical and Astronautical Engrg., Univ. Illinois, Urbana, 
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which occur different parts the building based largely the experience 
the analyst. This may lead errors, since damping has over-riding effect 
the structural response. 

the present paper, alternate formulation used which dispenses 
the intermediate step computing the normal modes and deals directly with 
both displacement and force types structural response. Making use the 
concept transfer matrix (8), the approach allows realistic representation 
the interaction between various parts within the system. particular, 
possible assign appropriate damping coefficients associated with various 
structural members joints between members, both. The approach especially 
suitable for tall buildings having many stories the construction. Additional 
simplification achieved under assumption that every story the building 
identically constructed, although this assumption not required for the validity 
transfer matrix approach. The application the methodology illustrated 
example 40-story building under the excitation wind turbulence 
field which statistically stationary time but nonhomogeneous space. 


The structural model chosen for the present study N-story building 
shown Fig. The following assumptions are made simplify the analysis: 


(c) 


FIG. 1.—Structural Model: (a) N-Story Building under Wind Loads; (b) Construction 
Unit; (c) External and Forces Unit 


(1) The Building structure composed identical construction units typified 
the one shown Fig. (2) the mass each unit concentrated 
the floor level; (3) linear elasticity provided massless columns suear 
walls between neighboring floors; (4) wind excitations are applied discrete 
floor levels; and (5) the response the structure external excitations can 
described displacement variable and force variable each construction 
unit. 

Referring Fig. let and Y,_, be, respectively, the displacement 


the floor and the shear force the column shear wall the unit. 
seen that 


(a) (b) 
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which the mass the jth floor; damping coefficient; and 
the stiffness. 


Taking the Fourier transform and rearranging terms, these two equations may 
cast matrix form follows: 


or, more simply: 


j 


which {Z}, {X,, state vector; over bar denotes the Fourier 
stiffness complex-valued the structural damping viscous damping 
between floors are taken into consideration. The physical meaning matrix 
[T] clear; represents the transfer mechanism construction unit. 
the excitation absent, state vector transfered state vector 
through 


Eq. can applied repeatedly obtain relation between {Z}, and 
namely: 


The boundary conditions the base and the top floor the building dictate 
that 


Substituting Eq. into Eq. obtain from the second row 
N 
j=l 
element [7] follows from Eq. that 


Thus the vector (Eq. expressible terms the random excitations 
and the elements transfer matrix [7]. 

The Fourier transforms the displacement the mth floor and the shear 
force immediately above the mth floor, X,, and which constitute the response 
vector {Z},,, can calculated the same manner Eq. With being 
replaced Eq. becomes 
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Substituting {Z}, {0, and using Eq. obtain 


Eq. input-output relation the (transformed) frequency domain 
which the ... are the inputs and the outputs. The various 
elements characterize the structural model. 


When wind loads building are modeled stochastic processes, the 
inputs Eq. are Fourier transforms stochastic processes. 
the wind loads are assumed statistically stationary time, then their 
Fourier transforms not exist the sense mean-square integrals. this 
case, the are merely formal representations; however, Eq. can still 
used conveniently construct spectral relationships between inputs and outputs. 
For example, the spectral density the random displacement X,, the mth 
floor, denoted can obtained multiplying the first row 
Eq. its complex “conjugate, and replacing each product F,F* the 
cross-spectral density namely: 


N N 


kel 


i=l kel 


which Re{ signifies the real part complex expression. Similarly, the 
spectral density the shear force can obtained from the second row 
Eq. 10: 


N N 


kel 


Cross-spectral densities the outputs such (w) can also obtained 
desired. 

seen from the preceding equations the computation output spectral 
and cross-spectral densities requires the knowledge and The 
former are the cross-spectral densities the wind loads, and the latter are 
the elements [7]”. These will considered the subsequent sections. 
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Analogous the case airfoil, both drag and lift forces are generated 
due relative motion between air particles and building. The drag forces, 
however, are responsible for the along wind motion the building. Since most 
buildings are not streamlined, the actual force-generation mechanism compli- 
cated and not fully understood the present time. For lack better 
representations, the common practice has been the use rather simple analytical 
form for the drag forces follows: 


which air density; tributary area the jth floor; wind 
velocity relative the structure; and drag coefficient which can different 
for different structural forms and should experimentally determined, least, 
for each category buildings. 

The velocity the building itself generally low. Thus can taken 
just the wind velocity, which may decomposed into mean velocity 
and random fluctuation about the mean. Let 


which assumed stationary random process with zero mean. 
Then E[V,] where denotes ensemble average. reasonable 
13. Thus 


has been found (1) that the mean wind velocity varies with the height 
from the ground according power law 


which the height the jth floor; reference height called the 
gradient height; the wind velocity the gradient height; and constant 
between 0.15 and 0.5. Roughly speaking, the gradient height one beyond 
which the mean wind velocity becomes nearly independent the distance from 
the ground. This parameter, well the exponent Eq. 17, related 
the ground roughness, characterized the so-called roughness length which 
better known the fluid dynamicists. Davenport (3) has provided chart 
where and can found from the roughness length. 

should noted that have identified the mean wind velocity being 
the same the ensemble average which constant for stationary random 
process. Several writers have taken, perhaps, broader point view treating 
the mean wind velocity random variable. particular, Davenport (5) has 
suggested that the mean wind nearly Rayleigh distributed, and that the extreme 
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value the mean wind given location follows the Type extreme value 
distribution. However, for structural response calculations, the two approaches 
lead the same results for structural response statistics. 

The autocorrelation function given 


which the autocorrelation function for terms spectral 
densities, Eq. equivalent 


which denotes the Dirac delta function. Following Davenport (4), the 
power spectrum the fluctuating part the wind, has nearly universal 
expression 


which frequency, hertz (cycles per second); reference mean 
wind velocity, meters per second above the ground; constant; 
and 1,200 n/u,, cycles. 

Equation one-sided spectrum the positive frequency domain. When 
converted correspond the nondimensional and two-sided spectral 
density Eq. becomes 


interest note that the Davenport spectrum tends zero 
for large similar the well-known von Karman spectrum, but they differ 
the low-frequency region. The von Karman spectrum becomes rather flat 
zero-power law) diminishes, whereas the Davenport spectrum tends 
linearly zero with Furthermore, the Davenport spectrum implies that the 
scale turbulence independent height, similarly proposed Berman 
(2), although intuition would suggest that the turbulence scale should grow with 
height. The preceding brief description bring out subtle implications 
Eq. which has been used our numerical calculations. The methodology 
presented this paper not restricted any specific input spectra. 
not our intention make exhaustive comparison suggested spectra 
the literature, some which have not been cited herein. 

Analogous and 19, the cross-correlation function and the cross-spec- 
tral density and are, respectively: 


Following Davenport (4), the one-sided cross-power spectrum for and 


defined the positive domain (in cycles per second) may expressed 


S(n) exp 
u, 


/ S(n) 4k,u? x? (20) 
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story height. Converting Eq. for the nondimensional and and for 
two-sided cross-spectral density obtain 


Equation reduces Eq. when Equations 16, 17, 19, 21, 23, and 
provide the required statistical information the wind loads for computing 
the first and the second statistical properties the structural response. 


appearing Eqs. 11, 12, etc. for the spectral and cross-spectral 
Their computation the obvious chain multiplication the matrix 
times can become very tedious when large, and the result can 
inaccurate due accumulated roundoff errors. such case, advantageous 
use the following alternate computational scheme which based some 
special properties common all transfer matrices (6). 

Recall that 


eer 
can easily shown that |7| and that the characteristic equation for 
either the eigenvalues the eigenvalues [7] the same 


Therefore, the roots Eq. are reciprocal pair which may expressed 


or, denoting K/m and B/(2mw,): 


known that any analytic function matrix may expressed 
linear combination any linearly independent analytic functions the 
same matrix. choosing the form 


n a, b, 
the coefficients and can determined using the fact that Eq. holds 
when the matrix [7] replaced any eigenvalue [7]. The two relationships 
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resulting from substituting [7] exp lead 
cos 
cos 
which, turn, lead 


The elements and appearing Eqs. and are given specifically: 


sin 


should noted that the value for complex for damped system, 
when the system undamped where can seen from 
Eq. 29. these cases, can computed from 


which the right hand side Eq. 29. 

Numerical Examples.—To illustrate the application the present methodology, 
and are used Eqs. and compute the spectral densities 
for displacement response and shear force response. each case, the result 
can separated into two parts; i.e.: 


which 


N 


j=l 


w h m m 


u, k=l 


r i=l kel 


‘ 
sin 
(38) 


EM2 ALONG-WIND MOTION 


N m 


sj 


clear from and that and n,, are the average displacement 
and average shear force 


whereas and gives the frequency distribution the variances; 


40-story building model has been chosen for the numerical example. The 
structural and aerodynamic physical data used the computation are follows: 
Number stories 40; individual story height lumped mass 
individual floor 1,290,000 kg; elastic stiffness each story 10° 
N/m; damping coefficient 21,550 N/m/sec; wind-load tributary area for 
each story 192 m’; gradient height 300 gradient wind velocity 
44.69 m/sec; reference mean wind velocity height 11.46 
m/sec; drag coefficient 1.2; air density 1.23 kg/m’; ground surface 
drag coefficient 0.03; exponent for the mean-wind-profile power law 
0.4; and 7.7. 

Figure shows the computed spectral density for the displacement the 
top floor, without the dirac delta function Therefore, this figure 
actually graph the frequency distribution the variance, i.e., the 
function given Eq. for the special case 40. expected 
that the greatest displacement excursion occurs the top floor. the other 
hand, the maximum shear force occurs most likely the bottom story. Thus, 
illustrate the nature the shear force response, the frequency distribution 
the variance plotted Fig. This the function given 
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Eq. 40, computed for the special case and 40. Only the positive 
domains are shown Figs. and the graphs for the negative domains 
can obtained their mirror images. Two peaks are seen both figures. 
The first peak 0.04 rad/sec due the maximum the wind spectrum, 
whereas the second peak 1.02 rad/sec coincides with the fundamental natural 
frequency the structural model. This means that both and are dominated 
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FIG. 2.—Frequency Distribution Variance Displacement Top Floor 40-Story 
Building 
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FIG. 3.—Frequency Distribution Variance Base Shear Force 40-Story Building 


the fundamental mode; however, higher modes may have greater contribution 
toward the displacement force response different locations the structure. 
Variances and were obtained integrating numerically their 
frequency distributions (Eqs. and 43). Upon taking the squareroots, the standard 
can compared with the respective mean values 0.1013 and 
0.389 obtained, respectively, from Eqs. 37, 39, and 41. 


1 
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number simplifications have been used this initial attempt analyzing 
the wind-induced vibrations multi-story building with transfer matrix 
formulation. restricting the building motion the along wind direction and 
considering only one displacement variable and one force variable each 
story unit, only transfer matrices are required represent the structural 
model, permitting closed form solutions for some the results. The order 
the transfer matrices will increased when the torsional and the cross-wind 
motions are taken into account, subject which will dealt with paper 
appear later. 

The assumption that every story the building identically constructed 
has resulted substantial saving computer time. Specifically, long chain 
product transfer matrices replaced the basic transfer matrix raised 
high power, and simple scheme has been devised facilitate the 
computation. The use essentially identical construction for every story 
actually rather common practice the building industry because the overall 
cost effectiveness, least for moderately tall buildings. For very tall buildings, 
may more economical divide the whole building into several groups, 
each which consists array identically constructed units. the latter 
case, the method described herein for computing arbitrary power transfer 
matrix can still used advantage. 
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The following symbols are used this paper: 


tributary area floor; 

drag coefficient; 

constant appearing Eq. 24; 

wind force applied jth floor; 

story height; 

stiffness one story; 

constant appearing Eq. 20; 

mass jth floor; 

average value X,,, Eq. 41; 

average value Eq. 41; 

frequency distribution variance 
frequency distribution variance 
autocorrelation function 
autocorrelation function 

real part complex expression; 
one-sided spectrum fluctuating part wind velocity; 
transfer matrix; 

reference mean wind velocity meter above ground; 
mean wind velocity floor; 

wind velocity gradient height; 

relative wind velocity jth floor; 
displacement floor; 

shear force above floor; 

height jth floor; 

reference height; 

jth state vector; 

height separation /)h; 

damping coefficient; 

frequency, radians per second; 

air density; 

constant appearing Eq. 17; 

dirac delta function; 

mean wind force applied floor; 
eigen value [7]; 


Cp 
F, = 
‘e{ 
u, 
xX, 


ALONG-WIND MOTION 


exponent 

random fluctuation wind velocity about mean; 

spectral density 

spectral density F,; and 

over-bar quantity, indicates Fourier transform. 
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BUILDING INTERNAL PRESSURE: SUDDEN CHANGE 


The safety buildings under the action strong wind often depends 
much the building internal pressure the external pressure. Internal 
pressure also affects air infiltration and heat loss through buildings, indoor air 
pollution, occupant comfort, spread fire and smoke buildings, rain penetration 
into buildings, proper operation certain mechanical equipment buildings, 
and the accuracy indoor barometric measurements. Despite all these implica- 
tions, internal pressure has received far less attention than external pressure 
researchers, and consequently, the subject usually poorly understood. 
This paper seeks clarify the nature internal pressure variations, and 
provide rigorous theory the subject based fluid mechanics. 

When aperture building, such window door, opened suddenly 
wind, rapid change pressure happens inside the building, causing the 
internal pressure reach equilibrium steady-state value within short 
period. The equilibrium internal pressure may predicted the manner 
considered the first writer (4), and Newberry and Eaton (7). The transient 
response the internal pressure which happens before the pressure reaches 
equilibrium has been considered briefly Euteneuer (1) and the first writer 
(5). the intent the present paper provide more complete treatment 
and analysis this problem. related problem, namely, the fluctuations 
building internal pressure induced wind turbulence, has been treated recently 
Holmes (2). 

believed that correct prediction the transient response the internal 
pressure essential many facets wind engineering, such the analysis 
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the dynamic windload buildings—especially claddings. Questions such 
whether gust factor should used building codes determine internal 
pressure also can answered more intelligently when this transient response 
problem clearly understood. 


TREATMENTS 
The problem treated the transient response the pressure inside 
building following the breakage window door during wind storm. 


The window door may broken either the inequilibrium the internal 
and external pressures, wind-generated missiles, any other cause. 


(a) BEFORE BREAKAGE (t< 


(b) AFTER WINDOW BREAKAGE 


FIG. 1.—Change Pattern and Internal Pressure Due Window Breakage 


Referring Fig. suppose window the windward side broken 
time and suppose the wind speed The external and internal pressures 
are respective and p,,. soon the window broken, the 
internal pressure will rise, whereas the external pressure will remain constant. 
will take short while before the internal pressure will reach its equilibrium 
value which the same the external pressure p,. The problem can solved 
assuming either isothermal adiabatic condition, with the latter being 
more realistic. 

Isothermal Quasi-Steady Solution.—Perhaps the simplest solution the pre- 


viously mentioned problem use the Bernoulli equation for steady incom- 
pressible flow: 


BUILDING INTERVAL PRESSURE 


which velocity; pressure; and external air density. Subscripts 
and p,. Also note that for simplicity, the density the air inside the 
building has been assumed the same outside. Since much smaller 
than V,, and since the volumetric discharge where the effective 
area the opening 


2 


From the conservation mass 


which the total mass the air the building; the volume 
the building; and the density the internal air. 


Eliminating from Eqs. and using the equation state perfect 
gas, and assuming constant temperature 


which the ambient pressure; the stagnation pressure 


From the equilibrium condition reached when 


t, eq Pa P. 


Isentropic Quasi-Steady Solution.—A more refined approach assume 
isentropic, i.e., adiabatic and frictionless, flow and give due regards the 
variation density from 1-2. The Bernoulli equation for this case (6): 


which the ratio the specific heat capacity constant pressure and 
constant volume, i.e., C,/C,. The continuity equation 


dp, 
(7) 


The equation adiabatic gas 


P2 Pa 


which the subscript refers the condition. 
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Substituting Eq. into Eq. and rearranging terms yields 


V k 1/k 


Using the transformation and t/t, Eq. changed to: 
Integrating Eq. and using the initial condition when yields 


(k—1/k) 


The Eq. represents natural logarithm. 

Isentropic Unsteady Solution.—The aforementioned isentropic solution the 
quasi-steady solution because the Bernoulli equation used, i.e., Eq. for 
steady flow. get the true unsteady isentropic solution, unsteady term 
must added Eq. which now becomes (6) 


which infinitesimal distance along streamline. Note that the term 


has been dropped from Eq. because V,. 
The line integral Eq. cannot evaluated exactly. However, approximately 


which the effective length the jet entering the building. 

The previous simplification possible because, shown Fig. 1(b), the 
air jet entering the window opening contracts only slightly. forms jet column 
which can approximated cylinder having length and cross-sectional 
area The length usually the neighborhood 0.8d, which 
the diameter the opening (if circular), the length each side the 
opening (if squared). 


yields 
(10) 
Ps; 
dV, 
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Substituting Eq. into Eq. and using Eq. can shown that 


Next, using Eq. and the conditions given Eq. 11, the previous equation 
becomes 
2AL, 2AL, 
Ordinarily, V,/2AL, and the aforementioned equation reduces 
2AL, 


The initial conditions for solving Eq. are 


Comparison 


Results the three foregoing solutions can compared example. 
The example uses the following values: 1.1 m/s; 
259 (from Eq. 20). 

From the foregoing values, Eq. becomes 


t, 


becomes 


and become, respectively 
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The three different solutions are plotted Fig. The solution Eq. 
was obtained numerically using the standard Runge-Kutta method with different 
integration intervals. was done using CSMP (Continuous System Modeling 
Program) Amdahl 470/V7 computer. Note that small but noticeable 
differences exist among the three solutions. 

can seen that initially the rise the internal pressure predicted 
Eq. not fast the rise based the quasi-steady isentropic solution. 
This due the inertia effect the flow not considered the quasi-steady 
approach. can also seen that the internal pressure becomes the same 
the external pressure time roughly equal the time constant 
t,. Since the time constant this example sec, the internal pressure 


2 
— - 2 - o4 s 


%,70 155 LOG (: 165 


(ISENTROPIC) 


4 


0.010 


FIG. 2.—Transient Response Internal Pressure 


the building reaches the external pressure within 0.2 sec the window 
breakage—a rather fast response indeed! 


INTERNAL 


The preceding study the change the internal pressure following window 
failure was carried out until has reached the value was assumed 
implicitly that equilibrium, i.e., steady state, has been reached soon 
has risen p,. reality, the inertia air will cause exceed 
before the flow stopped. Then, the higher pressure inside the building causes 
the flow reverse direction which turn causes decrease air moves 
out the building. The reverse flow does not stop until becomes somewhat 
less than p,. After this, air rushes back into the building the beginning 
the first cycle. The flow will oscillate this manner until either the fluctuations 
are damped out after many cycles oscillation building failure occurs. The 
phenomenon similar the water hammer waves caused the sudden closure 
texts fluid mechanics. 
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The aforementioned oscillatory motion has been treated Holmes (2) 
1979 using the classical Helmholtz resonator model studied acoustics (8). 
The Helmholtz model based assumed analogy between the acoustic 
resonator and mechanical system consisting mass attached spring 
involved damped oscillations. Based the Helmholtz model, Holmes derived 
the following equation predict internal pressure variation (2): 


which the polytropic exponent which equal 1.4 and unity, respectively 
for adiabatic and isothermal flows; the orifice coefficient; the equivalent 
length the slug air the opening; the opening area; and dot and 
double dots above any quality represent respectively the first and second 
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Vv =30 m/sec 
Vo 
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Pio *Pamb *!01300Pa 
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HOLMES’ EQUATION 
(c=06) HOLMES EQUATION 
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FIG. Variation Internal Pressure Due Sudden Window Breakage 


derivatives with respect time, The four terms Eq. are, consecutive 
order, the mass-inertia, viscous-damping, elastic-stiffness, and forcing-function 
terms. 

will now shown that equation similar Eq. can derived 
from the Bernoulli equation using the same approach presented before 
the section Isentropic Unsteady Solution, except for the fact that reverse 


flow assumed whenever negative. Using such approach, the Bernoulli 
equation yields 


which the jet velocity which the same for forward flow. 
The plus and minus signs the last term Eq. correspond 


a“ 
7 1 Vi \ 
v 
° 
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(reverse flow) and (forward flow), respectively. 
Using approach similar the derivation Eq. 18, can proved 
that Eq. leads 


i 
2kAL, 
Knowing that the sign the last term 


Ao =1.0 

m/sec 

Vo =600 m3 

Pe 

Pio *Pamb =!01300pe 
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FIG. 4.—Isentropic Variation Internal Pressure Due Sudden Window Breakage 


p,, i.e., plus when negative and on, the term may written 


Substituting Eq. into Eq. 28a yields 


z z 1 E 


Equation the counterpart Eq. Holmes. Both equations may 
used for the determination the variation the internal pressure caused 
external pressure change. Note that both equations are second-order nonlinear 
ordinary differential equations. Since Eq. has more rigorous foundation, 
believed more correct than the Holmes equation. can shown 
that using the simplifying assumption p,, the same derivation 
which led Eq. yields Eq. 26. Thus, Eq. special case Eq. 
when the density change pressure change small. 


k-1 
and 
28a always opposite the sign 
HOLMES’ EQUATION 
(e=06) 
| / 
0.50 150 2.00 
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numerical example will now considered compare results computed 
from that from the Holmes equation. Consider the change internal 
pressure following sudden window breakage. Assume the external pressure 
the window stagnation pressure 1), and the initial internal pressure 
the same the ambient pressure which atmospheric. Furthermore, assume 
wird speed, building volume, and the window area respectively m/s, 
600 and These numbers are the same used Holmes (2) 
example. contraction coefficient 0.61 used herein obtain the effective 
area from the window area A,, i.e., 0.61 A,. Both and are assumed 
this case. Values C,, versus time were computed using Eq. and the Holmes 
equation, and the results are plotted Fig. Although Holmes used value 
1.2 (between isothermal and isentropic), value 1.4 was used this 
analysis, the flow was assumed isentropic. The results Holmes’ equation 


HOLMES 


/ / \ 


(SEC) 


FIG. 5.—Effect Variation Internal Pressure Based Eq. 


this analysis not differ much from those presented his paper using 
1.2. 

shown Fig. C,, versus based Eq. similar that based 
Holmes’ equation with orifice coefficient 0.6. For this case with 
large window, although the equilibrium internal pressure coefficient the same 
the external pressure coefficient which 1.0, the maximum C,, during the 
first cycle oscillation about 1.8. This means overshooting internal 
pressure much 80%. 

one uses the same data the foregoing example, except window 
area 1.0 which much smaller than the one used the example, 
the same computations lead the results given Fig. Note that the maximum 
overshooting for this case with smaller window less than 25%, and the 
damping faster than the previous case with large windows. 

interest see what value agrees the best with the Holmes result based 
0.8d. For this purpose, C,, versus was computed based values 


/ / / \ / 
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ranging from The results are plotted Figs. and respectively 
for large and small windows. Note that both cases the computation based 
Eq. agrees closely with Holmes’ results when approximately 0.5d. 
The frequency period oscillation both cases changes with smaller 
gives higher frequency smaller period. 

The period oscillation calculated from Holmes’ equation with 0.8d 
0.23 sec and 0.40 sec, respectively for and These periods 
are approximately the same those computed from Eq. using 


SEC 


EQUATION 


SEC 


(SEC) 


0.5d. The corresponding periods the undamped flow Helmholtz resonator 
can calculated from (8) 


which was assumed that 0.61 A,. Using Eq. 31, the values 
for and are respectively 0.26 sec and 0.46 sec. 


ano 


This paper demonstrates that the variation the internal pressure building 
following sudden window breakage can predicted from the unsteady isentropic 
Bernoulli equation. The result second-order nonlinear differential equation 
which better approximation than the equation derived from the Helmholtz 
resonator given classical acoustics. While the result based the Bernoulli 
equation valid for small well large pressure variations, the result based 
the Helmholtz resonator model correct only for small pressure 

Two matters need clarifications. One that may appear that the Bernoulli 
equation was not valid for this case because frictional loss. This not 
true, however, because the Bernoulli equation was used only between point 
and point Fig. which are within region convergent flow with 
little energy loss. The bulk the energy loss between point and point 


Letl.od 
1.0 
Cpi 
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where the Bernoulli equation was not used. This loss was automatically taken 
means complete loss the kinetic energy the flow. Another matter 
concern the validity the implicit assumption that the building rigid. 
reality, building walls and roofs deflect somewhat under the action high 
pressure. This building flexibility may provide extra damping the flow. How 
severe this effect under ordinary condition not known. area that 
deserves attention and future research. Vickery (11) has suggested empirical 
approach which modified bulk modulus air used embody the building 
flexibility effect. 

The study concludes the following: 


Following the sudden breakage window door strong wind 
the order m/s, the internal pressure adjusts the external pressure 
level time approximately equal 0.02 which usually much less than 
second. The previously mentioned finding implies that when window 
door open during high wind, the internal pressure will fluctuate with the 
external pressure even when the gust frequency higher than Hz, unless 
the window unusually small, the building unusually large. The aforemen- 
tioned findings are general agreement with evidences reported recently 
other investigators 

After the internal pressure has reached the external pressure, the internal 
pressure does not cease rise. continues increase until the internal pressure 
becomes much higher than the external pressure. Then, reversal flow direction 
occurs, and the internal pressure starts fall. Large overshooting the internal 
pressure may produced when large window broken during windstorm. 
This may have serious implications building safety. 

Sudden breakage window causes damped oscillation internal pressure. 
The frequency this damped oscillation approximately the same that 
predicted for undamped Helmholtz oscillator. The frequency increases when 
the area the opening increases. 

Experimental data are needed determine the value before Eq. 
can expected yield accurate results. Likewise, experimental data are 
needed determine the orifice discharge coefficient before Holmes’ equation 
can expected yield accurate results. the absence experimental data, 
recommended that value between 0.5 and 0.8 used for computing 
internal pressure variations. 
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The following symbols are used this paper: 


effective area opening; 
actual area opening; 
defined Eq. 11; 

external pressure coefficient; 
internal pressure coefficient; 
initial value 
orifice discharge coefficient; 

specific heat constant pressure; 
specific heat constant volume; 
defined Eq. 

dimension square window; 

defined Eq. 285; 

distance along streamline; 

length air jet entering building; 
length air slug acoustic resonator; 
mass; 

polytropic exponent; 

pressure; 

atmospheric pressure; 

external pressure; 

internal pressure; 

initial value p,; 

stagnation pressure (p, 
volumetric discharge; 


Soak 


=< 
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time; 
time constant 
wind velocity; 
velocity air jet entering building; 
internal volume building; 
initial value 
initial value 
1)/k; 
ambient-air density; 
external-air density; 
internal-air density; 
period damped oscillation computed from Eq. 30; and 
period undamped oscillation calculated from Eq. 31. 
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TORSIONAL VIBRATION ALONG-WIND 
STRUCTURES 


INTRODUCTION 


Although the along-wind vibration structures subjected wind loadings 
has been active research topic for the past twenty years, very few investigators 
have considered the rotational vibration these structures. Torsional motions, 
not considered design, could result overstressing the columns 
the periphery the structure and could also damage cladding attachments. 
The majority the research develop methods analysis has been directed 
toward predicting the expected maximum along-wind translational response 
structures due buffeting the wind. The methods that have been developed 
range from equivalent static load methods (1,14) which utilize design charts 
more precise methods requiring computer analysis (12). 

Hart (7) derived general set equations that could used analyze 
the three-dimensional response buildings stochastic wind forces. The required 
experimental data was not available the time, however, apply them 
practical problems. Recent work Patrickson and Friedmann (10) indicated 
that the dynamic torsional response buildings may lead displacements 
that are the same order the along-wind motions. Their work also indicates 
that velocities and accelerations due torsional motion may significant. 
experimental and analytical investigation Hart, al. (8), demonstrated 
the importance torsional vibrations buildings subjected wind and 
earthquake loadings. Both the first and second writers (6) demonstrated analyti- 
cally that even perfectly symmetric structures may experience serious torsional 
motion when buffeted the wind. Several investigators not cited here have 
provided experimental evidence relative magnitude the wind-excited 
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torsional response symmetric well nonsymmetric buildings. more 
thorough historical review the subject may found Ref. 11. 

this paper, methodology for analyzing the torsional vibration wind 
excited structures presented. The method based random vibration concepts 
and yields the expected maximum translational and torsional responses. The 
responses several single-mass structures are investigated. The goals this 
research were twofold. The first objective was develop model for analyzing 
the torsional response structures that paralleled those used for analyzing 
along-wind responses. Thus the approach would have the advantage both 
being firmly founded the existing body knowledge regarding wind structure 
and its effects, and being easily adopted engineers familiar with current 
methods. The second objective was investigate the response several simple 
structures order identify which wind and structural properties significantly 
influence the torsional response and determine whether the predicted torsional 


A: Edge Point 
E: Elastic Center 
G: Moss Center 


FIG. 1.—Schematic View Problem FIG. 2.—Geometry Problem 


responses were large enough warrant extending the method building-type 
structures for use design. The results clearly indicate that torsional vibration 
can contribute significant amount the total motion wind-excited structure. 


Consider the schematic idealized single-mass structure with approaching 
wind shown Fig. For the coordinate system located the mass center 


and for the notation shown Fig. the equations motion the system 
may written 
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which F(t) and respectively, the random fluctuating force and torque 
that act the structure; and and respectively, the structure’s translational 
and rotational stiffness with respect the elastic center. general, the equations 
motion are coupled. addition, will shown that F(t) and may 
correlated. Since F(t) and are random variables space and time, 
the dynamic responses y(t) and will also random variables. 

the frequency domain the force-response relationships may written 


which the system transfer function; and and [S,(n)] 
respectively, the spectral density matrices the random forces and the system 
responses. 

The derivation the system transfer function (sometimes called the mechanical 
admittance frequency response function) well known and may written 


H = 


which [M], [C], and [K] respectively, the system mass, damping, and 
stiffness matrices; and Since all the elements [H(n)] are readily 
determined from the properties the structure, one needs only derive the 
elements the input spectral density matrix able determine the spectral 
density matrix the response. Although the relationships for the translational 
vibrations are well known, they will included the following derivations 
for completeness. 


Wino Forces Structures 


The wind velocity, may assumed composed the sum 
its mean and fluctuating parts: 


which v(x, considered zero mean, normally distributed stationary 
random process. Neglecting terms that contain derivatives the wind velocity, 


further assumed that the pressure point the structure’s face 
given 


which the density air; and pressure coefficient. Since the 
fluctuating part the velocity small compared its mean, the fluctuating 
pressure the face the structure may approximated 


These fluctuating pressures produce the randomly varying force and torque 


A 


(5) 
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which the sum the averaged pressure coefficients the windward 
and leeward faces. The use this manner implies that the pressure 
fluctuations the windward and leeward faces are perfectly correlated. Even 
though experiments have suggested that this assumption not usually true 
(9), results conservative estimates displacements (12), greatly simplifies 
the following derivations, and probably not unwarranted light all 
the other assumptions that are made. also implied Ref. that pressures 
the across-wind faces produce negligible torque. This would normally 
true the along-wind dimension the structure small compared its width 
and the structure’s stiffness and mass are symmetric about its axis. 

The entries the spectral density matrix the forces are 
and the cross spectra, which are respectively the Fourier transforms 
the correlation functions the force and torque and their cross-correlation 
function. These may written terms the cross spectra the fluctuating 
velocity 


4 a 


A A 


4 4 


Based experimental evidence, the cross spectra the wind velocity may 
written the product the spectrum the wind velocity and the coherence 
function (14) 


which and the exponential decay coefficients for the and directions 
respectively. Using Davenport’s velocity spectrum (2) 


Simiu (12) proposed equation for the wind spectrum that believed 
better founded theory. However, Eqs. and have been more commonly 
used. should pointed out that this paper concerned primarily with the 
relative magnitudes the along-wind and torsional vibrations. Consequently, 
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the influence the various assumptions regarding the velocity spectrum and 
the pressure correlations the results should minimal since their effects 
should roughly equal for both types response. thorough analysis 
most the assumptions and approximations used analyzing the response 
buildings wind loadings may found Ref. 13. 

Thus far, the development the random wind forces structure has 
been general and would apply any structure that rectangular plan. The 
computations may greatly simplified for the simple two-degree-of-freedom 
structure Fig. Since the mean wind velocity varies slowly with height 


which the height the center the structure. Thus the mean static 
force the structure 


which the width; and the depth the structure. Using the relationships 
presented, Eqs. 10-12 may written 


4F? 
2 2 


2 


which and the nondimensional aerodynamic admit- 
tance functions which are derived follows. The spectral density matrix 
the excitation now defined and Eq. becomes 


ey Soo Ho 


WA, 


Functions 


nondimensional variables and z/D are defined, the explicit 
forms the aerodynamic admittance functions used Refs. 18-20 become 


7 
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and 


1/2 


general, the evaluation the admittance functions requires four-fold 
numerical integration. However, these may evaluated analytically use 
the approximation suggested Davenport (4): 


Performing the integrations yields the following analytical expressions for 
the aerodynamic admittance functions: 


ne 3D 1 2 2 1 2 4 


x 


x z 


The aerodynamic admittance function for translation, Eq. 30, was developed 
previously Davenport (1). These functions make possible the direct computa- 
tion the complete spectral density matrix the structure’s response. will 
demonstrated that this required for accurate estimate the expected 
maximum structural deformations. These equations will considered follows 
with regard specific example. 
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Once the functions the spectral density matrix the response have been 
formulated, the mean-square translational and rotational motions well their 
correlation coefficient may computed 


and 


Note that since F(t) and T(t) have been defined zero-mean random forces, 
the responses and also have zero mean values. 


assumed that the total displacement point the structure, 
independent height, the translation any point may written 


which Y,(x) the displacement the point due the static mean wind 
force, i.e.: 


which the mean displacement the coordinate center; and 
mean rotation. The rms displacement any point 


which the spectral density function the total response 


The expected maximum value time interval (usually assumed 
one hour) can written 


which peak factor developed Davenport (3), and approximately 
equal 


0 
0.577 
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Thus, estimate the most probable maximum displacement including torsion 
any point the structure can made. 


The equations derived were used analyze the responses several variations 
the structure Fig. The basic structure has width, equal 
400 (37.2 The reference wind velocity miles/h (129 km/h) 
and the translational natural frequency, n,, 1.0 were held constant for 
all examples. The effects magnitude the torsional natural frequency, the 
structure’s width, and the location the centers mass and rigidity the 


mph 
We s 


W=30ft 


— 
xix 


White Noise Velocity Spectrum 


FIG. 3.—Rotational Vibration Symmetric Structures 


maximum probable response values were examined. For each example where 
the structure nonsymmetric, the expected maximum edge displacement, max 
Y,, computed using Eq. normalized the expected maximum displacement 
that would computed the dynamic torsional motion were ignored. The 
latter quantity simply the expected maximum along-wind displacement for 
the center symmetric structure plus the displacement due the static 
rotation, The normalization was done this way because designers would 
normally compute the mean rotation their calculations. Thus, the results 
reflect the consequence ignoring only the dynamic torsional response design. 

The examples may divided into three classes according the form 
21. For perfectly symmetric structure the centers mass and resistance 
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both lie the vertical center line the structure’s face. For this case all 
the off-diagonal terms Eq. are identically equal zero. Thus, the 
equations motion are uncoupled. This means that the translational and rotational 
motions are statistically independent and thus may evaluated independently. 


Note that torsional vibrations will still excited, however, due the spatial 
randomness the wind. 


Frequency, cps 
(a) 


Wi=i6, 


Frequency, cps 
(b) 


Frequency, cps 


FIG. 4.—Aerodynamic Admittance Functions: (a) Translation; (b) Torsion; and (c) 
Cross 


The effects varying the structure’s torsional natural frequency the 
expected maximum displacement point are shown Fig. decreases 
with respect n,, the motions the edge the structure due torsional 
vibrations become large and may actually substantially larger than those 
due translation. Most the effect associated with lowering the torsional 
frequency results from the shape the wind spectrum which increases with 


0.10 
206 
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decreasing frequency this region the spectrum. This also demonstrated 
which was assumed excited artificial wind with ‘‘white 
spectrum, but with the same spatial correlation the natural wind. 

should also expected that increasing the width the structure would 
increase the rotational response. This should result because pressures the 
extreme edges the structure would less correlated and would lie further 
from the elastic center. Both these effects would produce greater dynamic 
torque. This demonstrated Fig. for structure with width and depth 
equal 30.0 (9.15 and 13.33 (4.07 respectively. These values were 
chosen that the total exposure area would remain constant. 

The second class structures has only one type asymmetry, either geometric 
structural. Geometric asymmetry occurs structure has its centers 


wW=D=20ft 


max X 4 (incl. rot. vib.) 
max (no rot. vib) 


2 


FIG. 5.—Correlation Coefficient Re- FIG. Vibration Struc- 


sponses Structures with only Geo- tures with only Geometric Asymmetry 
metric Asymmetry 


mass and rigidity the same line, but not the center line the exposed 
face. The equations motion for this case are still uncoupled, e.g., 
H,, However, the cross aerodynamic admittance functions, 
are longer zero. Therefore, the correlation coefficient for the responses 
also nonzero and coupling the responses results. The aerodynamic admittance 
the coordinate system located ft, ft, and (3.05 3.66 
4.27 and 4.88 from point are shown Fig. Note that invariant 
with respect the coordinate center and that and increase the 
geometric offset increases. Thus, the correlation coefficient coordinate-system 
dependent. Also, was previously mentioned, identically equal 
zero for (3.05 which corresponds the coordinate system 


40, 
lops 

os * 

4 


EM2 TORSIONAL VIBRATION 333 


located the geometric center the structure. The correlation coefficients 
for these same structures with varying torsional frequency are shown Fig. 
The highest correlation the responses occurs when should 
expected. This value approaches 1.0 the geometric offset increases. 
general, the response correlation increases approaches and W,/W 
increases. 

The effects varying the rotational natural frequency and W,/W the 
maximum displacement point are shown Fig. The results show the 
same trend for the perfectly symmetric structures. The torsional natural 
frequency has the largest effect the rotational response. addition, the 
rotational response increases W,/W increases. Note that for structure 
with its centers mass and resistance offset only 10% its width from the 
geometric center 0.6), the dynamic part the torsional displacement 


40. 40 


xix 


folie) Q20 025 
W/W 


FIG. 7.—Rotational Vibration Struc- FIG. 8.—Rotational Vibration Struc- 
tures with only Structural Asymmetry tures with both Structural and Geomet- 
ric Asymmetry 


the edge the structure not computed the maximum displacement will 
underestimated about 25% for and about 80% for 
0.5 n,. might expected that the highest responses would result for 
since the force and torque are more highly correlated for this case. However, 
the lower correlation between the force and torque for structures with 
overshadowed the fact that the spectrum the torque higher 
the lower frequencies. 

the mass center located the structure’s geometric center but the elastic 
center not, structural nonsymmetry results and longer zero. 
The cross-aerodynamic admittance functions will again zero, however, since 
the reference system located the geometric center the structure. Results 
for structures with different torsional frequencies and eccentricities are shown 


30 30 
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Fig. These results are very similar those for structures with only geometric 
nonsymmetry shown Fig. Neglecting the dynamic part the torsional 
displacement for structure with eccentricity 10% its width will once 
more lead underestimates the maximum total displacement 25% for 

The most general case occurs when the centers mass, resistance, and 
geometry lie different lines. This results all the matrices Eq. 
being full. Shown Fig. are the results for structures with their mass center 
offset 10% the structure’s width [W, (3.66 m)] from its geometric 
center. expected, these motions increase with decreasing torsional natural 
frequency. Fig. somewhat misleading since shows decreasing torsional 
response with increasing eccentricity for less than 0.1. this region, 
the elastic center moves from the mass center for the geometric 
center for W,/W 0.1. The fact that the torsional motion decreases this 
region indicates the distance between the elastic center and the geometric center 
where the average pressure center lies more influential than the distance 
between the elastic and mass centers. For structure with (0.61 
the elastic center the geometric center. ignoring the dynamic torsional 
response for this case, the maximum displacement the edge the structure 
would underestimated about 25% for and 50% for 0.5 
n,. 
Summary 


method for estimating the expected maximum dynamic torsional response 
wind-excited two degree-of-freedom structure was presented. Aerodynamic 
admittance functions were derived which were used estimate the spectral 
density function the random torque and cross-spectral density function 
the force and torque acting the structure. These are required for estimating 
the structure’s translational and torsional mean square responses. Results for 
several examples indicated that, general, the dynamic torsional response 
increases the width the structure’s exposed face increases, the structural 
geometric eccentricity increases, and the torsional natural frequency 
decreases. was shown that for eccentricity only 10% the structure’s 
width, the total response structure can significantly underestimated 
the dynamic torsional response not included the analysis. The results 
Figs. 6-8 suggest that approximately linear relationship exists between 
the torsional dynamic response and the eccentricity. this true for the full 
three-dimensional problems, simple design rules could result. 

General design rules were not derived from these results since the structures 
studied were simple and limited number. The results indicate, however, 
that torsional vibrations may problem real structures. Consequently, 
the writers are extending the procedures developed this paper three-dimen- 
sional building-type structures and will include across-wind forces the analysis. 
addition the effects higher modes vibration and pressure correlations 
along the full height the structure, the influence assumptions regarding 
the shape the design wind spectrum and the form the along-wind and 
cross-wind pressure correlation functions the torsional response buildings 
will investigated. 


TORSIONAL VIBRATION 


This research was partially funded the National Science Foundation under 
grant NSF ENV 77-07190 and the University Illinois 
This support gratefully acknowledged. The reading the manuscript 


Pecknold and Wen also appreciated. The results this study 
were summarized Ref. 


Davenport, G., Loading Factors,’’ Journal the Structural Division, ASCE, 
Vol. 93, No. ST3, Proc. Paper 5255, June, 1967, pp. 11-34. 

Davenport, G., Spectrum Horizontal Gustiness Near the Ground High 
Quarterly Journal the Royal Meteorological Society, Vol. 87, London, 
England, 1961. 

Davenport, G., Distribution Largest Values Random Function with 
Application Gust Proceedings the Institute Civil Engineering, Vol. 
28, 1964, pp. 187-196. 

Davenport, G., Prediction the Response Structures Gusty Wind,”’ 
International Seminar the Safety Structures under Dynamic Vol. 
Norwegian Institute Technology, June, 1977, pp. 257-284. 

Foutch, A., and Safak, E., Vibration Wind-Excited 
Proceedings, 3rd ASCE/EMD Specialty Conference, Sept., 1979, pp. 255-258. 

Foutch, A., and Safak, E., Vibration Wind-Excited Symmetric 
Journal Wind Engineering and Industrial Aerodynamics, Vol. No. 
Mar., 1981, pp. 191-201. 

Hart, C., Dynamics Due Stochastic Wind Journal the 
Structural Division, ASCE, Vol. 96, No. ST3, Proc. Paper 7154, Mar., 1970, pp. 
535-550. 

Hart, C., DiJulio, M., and Lew, M., Response High-Rise 
Journal the Structural Division, ASCE, Vol. 101, No. ST2, Proc. Paper 
11126, Feb., 1975, pp. 397-416. 

Lam Put, R., Response Tall Building Random Wind 
Proceedings, International Conference Wind Effects Buildings and Structures, 
1971, pp. 4.1.11. 

Patrickson, P., and Friedmann, R., Study the Coupled Lateral and Torsional 
Response Tall Buildings Wind UCLA-ENG-76126, University 
California Los Angeles, School Engineering and Applied Science, Dec., 1976. 

Safak, E., and Foutch, D., ‘‘Vibration Buildings under Random Wind Loads,”’ 
Structural Research Series No. 480, Civil Engineering Studies, University Illinois, 
Urbana, May, 1980. 

Simiu, E., Spectra and Dynamic Alongwind Journal the Structural 
Division, ASCE, Vol. 100, No. ST9, Proc. Paper 10815, Sept., 1974, pp. 1897-1910. 

Simiu, E., and Marshall, D., Loadings and Modern Building 
presented the April 22-26, 1974, ASCE National Structural Engineering Meeting, 
held Cincinnati, Ohio (Preprint 2268). 

Vickery, J., the Reliability Gust Loading Civil Engineering 
Transactions, Apr., 1971, pp. 1-9. 


The following symbols are used this paper: 


drag coefficient; 
pressure coefficient; 


D,,D, 
F(t) 


n 

P(x, 


y(x, n) 


T(t) 


V(x, 
v(x, 


z 
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exponential decay coefficients; 

viscous damping values; 

depth structure; 

defined Eqs. and 34; 

ensemble average; 

mean wind force; 

fluctuating total wind force; 

peak factor defined Eq. 43; 

height center structure’s exposure area; 
system transfer function; 

mass moment inertia; 

translational, rotational, and cross-aerodynamic 
admittance functions; 

surface drag coefficient; 

translational and rotational stiffness; 

mass; 

frequency, cycles per second; 

total wind pressure; 

fluctuating wind pressure; 

defined Eq. 29; 

spectrum fluctuating force, torque, and cross 
spectrum; 

spectral density matrices force and response; 
spectrum and cross spectrum fluctuating wind 
velocity; 

spectrum total response; 

spectrum dynamic translation, rotation, and 
cross spectrum; 

averaging time interval; 

fluctuating torque; 

reference mean wind velocity height ft; 
mean wind velocity center structure; 

total wind velocity point 
fluctuating wind velocity point 
width structure; 

distances from mass center edges; 

distance between mass center and elastic center; 
defined Eq. 15; 

cross-wind direction coordinate axis; 

total displacement height 

total displacement center height; 

mean displacement; 

mean displacement center; 

maximum probable displacement; 

along-wind coordinate axis; 

along-wind dynamic displacement; 

vertical coordinate axis; 

x/W; 
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mean rotation; 
dynamic rotation; 
frequency, radians per second; 
density air; 
coherence function, defined Eq. 26; 
Pye correlation coefficient; 
root-mean-square response for translation 
rotation; 
root-mean-square total displacement; 
coefficient given Eq. 28; and 
coefficient given Eq. 44. 
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NATURAL FREQUENCIES CURVED GIRDERS 


INTRODUCTION 


Curved girders are used frequently highway bridge structures. With the 
application digital computers and the development numerical methods, 
the mathematical difficulties associated with curved beam elements have been 
largely overcome. 

Extensive work has been published concerning the behavior horizontally 
curved members, summarized McManus, al. (14). These investigations 
are primarily concerned with the analysis curved beams subjected static 
loadings. However, comparatively little work has been done regarding the dynamic 
response curved beams. The free vibration circular ring has been studied 
Love (13) and Timoshenko (24), among others (15). Culver (6) and Chaudhuri 
and Shore (5) studied the free vibration horizontally curved beams. Practically 
all these previous papers lack generality; some are limited in-plane vibration, 
while others are restricted doubly-symmetric sections out-of-plane vibration 
only. 

The purpose this paper derive general matrix formulation determine 
the natural frequencies spatial curved beams shown Fig. including 
the warping contribution. Also included are the rotary inertia effects with respect 
flexure and torsion, and the effects antisymmetry cross sections, which 
were found important. 


Equations 


The equations equilibrium terms linearized componental deformations 
for the free vibration spatially curved beams projected coordinate axes 


Prof. Civ. Engrg., Marquette Univ., 1515 West Wisconsin Ave., Milwaukee, 
Wisc. 53233. 

Engr., Boeing Commercial Airplane Co., Seattle, Wash. 

Note.—Discussion open until September 1981. extend the closing date one month, 
written request must filed with the Manager Technical and Professional Publications, 
ASCE. Manuscript was submitted for review for possible publication March 18, 1980. 
This paper part the Journal the Engineering Mechanics Division, Proceedings 
the American Society Civil Engineers, Vol. 107, No. EM2, April, 1981. 
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FIG. 1.—Thin-Walled Curved Beam Coordinates 


shown Fig. can written (25) 


x 


R? 


which axial deformation the direction the circular axis, and 
displacements the direction the principal and axes, respectively; 
rotation about the shear center (primes denote derivatives with respect 
the axis, and dots represent derivatives with respect time); and 
elastic and shear moduli, respectively; mass density; polar radius 
gyration with respect the shear center; cross-sectional area; and 
coordinates the shear center taken from the centroid the section; 
and principal moments inertia the section; torsional constant; 
warping constant; radius curvature the circular axis; and 
coordinates the cross section; normalized sectorial coordinate; 
and and make pointwise right-hand Cartesian coordinates. 

The right-hand side terms Eqs. through represent equivalent static 
loading according D’Alembert’s principle. Eqs. and are coupled 
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terms and and govern the behavior the beam axial deformation 
and bending about the axis, usually the minor axis. Hence, the minor axis 
flexural behavior affected the extensibility the member. Similarly, Eqs. 


and are coupled terms and and describe the problem flexure 
and torsion. 


Minimum 


The variation the minimum potential energy functional for the curved beam 
can obtained modifying the expression obtained for straight beam (9,11). 
The internal forces acting the curved beam elements are 


which axial force; and internal bending moments about 
the and axes, respectively; St. Venant’s torque; warping torque; 
and bimoment. 

Examining Eqs. and the expression for the corresponding internal forces 
the straight beam (25) reveals that deformations projected coordinate 
axes the curved beam are coupled, due only the initial curvature. Thus, 
replacing the straight beam terms the variation the minimum potential 
energy functional for the straight beam with the corresponding curved beam 


terms (Eqs. results expression energy variation for the free vibration 
curved beam which 


N.= EBA -4) (2a) 
R 
1” 
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which J,(0,1) called the conjunct concomitant, which the sum 
all integrated terms. For natural boundary conditions, the conjunct equal 
zero and the system said self-adjoint. 


general displacement field for the static behavior curved beam 
assumed (27) 


which the asterisk denotes transposition. Each node rotates about three axes, 


(6) 


FIG. Forces and Displacements 


undergoes displacements along these axes, and warps, for total seven degrees 
freedom each node, shown Fig. order avoid unnecessary 
complication, the general displacement field divided into and pertinent 
the displacements associated with the two sets coupled equations, 
and and Egs. and Since these equations are coupled, there must 
coupled boundary condition each set equations. Introducing two 
kinematic boundary conditions (5), and which are defined 


+4) 
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the subdisplacement fields can written 


4,= 


which [N,] and matrices shape functions given Yoo (27). 
The generalized displacements are chosen nodal displacements, shown 
Fig. Thus 


which the subscripts and denote the left end and right end, respectively, 
element. The negative signs Eq. are due the assumed positive 
directions the corresponding degrees freedom shown Fig. 


For the dynamic problem under consideration, generalized displacement 
field assumed 


which A,(z) assumed displacement modes, i.e., the exact solutions 
the coupled homogeneous differential equations governing the static behavior 


the curved beam problem, Eq. The unknown amplitudes vibration may 
defined 


which angular velocity radians per sec, and Kroneckar delta, 


and all other (for j). Hence, the general displacement 
field for the dynamic problem can 


which {A} column matrix representing the nodal degrees freedom 
both ends element, and [N] matrix shape functions defining 
the elastic behavior the structure. The shape functions are used for both 
the static and dynamic problems. Therefore, the static solutions are exact 
regardless the grid refinement the structural model; the dynamic solutions 
were found extremely fast converging upper bounds, can expected 
due the use the variational procedure. 


Introducing the displacement field Eq. 13, and noting that 
Eq. may rewritten 


which and [m,] element stiffness and mass matrices, respectively, 
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the element. The elements the stiffness matrix, K,, are obtained from 
the first bracketed term Eq. and can expressed 


i 


Ny ” Ny, 
i ” ” 
” Nz, ” Ny, ” N ” Nw 


(in which 6,7, 11, 12, 13, 14) 


The elements the mass matrix, can derived from the second bracketed 
term Eq. and are expressed (9) 


Applying Castigliano’s Theorem, Part Eq. becomes 


which {Q} generalized forcing vector. For eigenvalue problems defining 
the natural frequencies, the vector becomes the null vector. Assembling the 
element matrix equation motion, Eq. 17, for the entire structure results 
the final eigenvalue equation for the dynamic problem 


Examination Eq. reveals that the element the mass matrix [m,] 
contains terms which reflect the effects the shear center location when 
does not coincide with the centroid the section. These effects are not reflected 
ordinary consistent mass matrix formulation. Therefore, the assembled 
structural mass matrix [m], when obtained setting shear 
center coincident with the centroid), will identical that obtained consistent 
formulation. Eq. general eigenvalue problem. 

combined Sturm sequence and inverse iteration technique, described 
Gupta (10), used for the solution the eigenvalue problem. The method 
efficient and numerically stable algorithm for the accurate computation 
specified roots and associated vectors the general eigenproblem. 
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order demonstrate the versatility and reliability the proposed method 
and the computer program developed, several example problems will presented. 
Some the results will compared with results obtained from previously 
published sources. 

Convergence.—The shape functions which describe the static behavior the 
curved beam were also used describe the dynamic behavior the curved 
beam. Therefore, the dynamic solutions will necessarily approximate. Due 


TABLE 1.—Convergence Single Span Curved Beam Solution 


Neglecting Effect Including Effect 
Rotational Inertia Rotational Inertia 


Difference, Difference, 


asa 
Number elements percentage percentage 


(1) (3) (5) 


A NEGLECTING EFFECT OF 
ROTATIONAL INERTIA 


OINCLUDING EFFECT OF 
ROTATIONAL INERTIA 


(RADS/SEC) 


8 
NUMBER OF ELEMENTS 


FIG. 3.—Convergence Single Span Curved Beam Solution 


the utilization the variational procedure, the dynamic solutions were expected 
upper bounds. The results the computer analysis confirmed this. 

Consider the example single span horizontally curved wide flange beam 
length 1600 in. (4064 cm). The cross-sectional properties the beam are: 

Values used for the elastic modulus and shear modulus are 29000 ksi 
(200.1 and 11200 ksi (77.3 The beam considered 
have all end conditions fixed. The radius beam curvature 1528 in. 


4.2106 +0.542 4.2078 +0.475 
4.1895 +0.038 4.1897 +0.043 
4.1881 +0.005 4.1883 +0.010 
4.1879 
420 
4.19 
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(3881.1 cm), and the beam spans angle 60°. The results the analysis 
are shown Table and Fig. 

The percentages variation are calculated using the result the 16-element 
configuration the exact value for comparative purposes. The convergence 
evident; even the result the relatively crude two-element configuration 
varies only minimally from the exact value. addition, can seen that 
the effects rotary inertia are quite insignificant, and can disregarded for 
practical purposes. 

Comparative Studies.—Culver (6) determined the natural frequencies 
horizontally curved beam utilizing the assumptions that the beam prismatic 
and doubly symmetric, and considering only vibrations normal the plane 
curvature. 

Shore and Chaudhuri (20) studied the free vibration horizontally curved 
beams, considering vertical displacement and rotation about the z-axis. Expres- 


TABLE 2.—Correlation Results for Curved Beam Frequencies 


Yoo and 
Fehren- Shore and 
bach Culver (6) Chaudhuri (20) 
Spanning Radius radians radians Difference, radians Difference, 
Angle, Curvature, per per asa per asa 
degrees inches second second percentage second percentage 


(2) (3) (4) (6) 


sions were derived for determining the natural frequencies curved beams 
both with and without the effects shear deformation and rotary inertia being 
considered. The previous problem presented here determined that the effects 
rotary inertia are negligible; therefore the expression for natural frequency 
curved beams without consideration the effects shear deformation and 
rotary inertia will used basis for comparison with the results this 
study. 

The example considered here single span curved beam length 
201.68 in. (512.27 cm). The beam curvature varied that the beam spans 
variety angles ranging from 10° 90°. Cross-sectional properties the 
5.042 in. (12.81 cm). The beam considered have pinned conditions 
both ends. Table the natural frequencies determined this study 


1155.5 204.7 202.5 1.1 203.3 0.7 
577.8 190.2 184.3 3.1 186.3 
385.2 165.8 162.2 164.7 0.7 
288.9 141.7 140.0 142.8 0.8 
231.1 121.3 120.9 0.3 122.8 1.2 
192.6 103.9 103.8 0.1 105.2 1.3 
165.1 89.6 88.9 0.8 90.0 0.4 
144.4 76.8 76.0 1.0 76.8 0.0 
128.4 64.0 64.6 0.9 65.2 1.9 


CURVED GIRDERS 


are compared with those calculated according Culver and Shore and Chaudhuri. 
Excellent correlation observed. 

Shear Center Effects.—Several examples were analyzed determine the effect 
the shear center location the dynamic response the curved beam. 
The first example considered was the singly symmetric section shown Fig. 
The curved beam has length 240 in. (609.6 cm) and has fixed conditions 
its ends. The radius curvature 275.02 in. (698.55 cm) and the beam 
spans angle 50°. Cross-sectional properties are: 15.4 (99.4 


SECTION 
PROPERTIES : 


15.4 in® 
480.48 

ly = 63.37 in* 

ly = 810 in® 

Ky = 1.93 in# 


o, = 0 
ay = —3.88 in 


EFFECTIVE MOMENT 
ARM OF INERTIAL FORCES 


CROSS SECTION DEFORMED SHAPE 


FIG. 4.—Singly Symmetric Curved Beam Cross Section in. 2.54 cm) 


CENTER 
CURVATURE 


TO CENTER 
OF CURVATURE 


ay = 1.634 


(b) 


FIG. Curved Beam Cross Section in. 2.54 cm) 


When the shear center effects are included, the resultant natural frequency 
the beam found 261.41 rads/sec. the shear center effects are 
eliminated (by setting 0), the natural frequency the beam 
262.97 rads/sec, only slightly stiffer than-the beam with shear center effects 
included. 

can seen the deformed shape the beam shown Fig. that 
the x-coordinate the shear center will contribute most significantly the 
inertial forces beam vibration. Therefore, for the singly symmetric section 
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symmetric about its vertical axis, the effect shear center location the 
dynamic response the beam will relatively small. However, beams 
not symmetric about the vertical axis, shear center effect will found 
significant. 

Two variations unsymmetric channel section are shown Fig. The 
cross-sectional properties the section are: (180.6 cm’), 
824.25 in.* (34305 180.05 in.* (7493.7 3661.6 (983286 
9.33 in.* (388.3 and 5.426 in. (13.78 cm). The beam 
has length 240 in. (609.6 cm) and spans angle 50°. Moments 
inertia and shear center coordinates are taken with respect the principal 
axes the section. basis for comparison, the shear center effects can 
eliminated setting the shear center coordinates equal zero. The results 
the computer analysis are shown Table 

Considerable variation seen due the effects the shear center, caused 
primarily the contribution the shear center x-coordinate. However, the 
effect the shear center location the dynamic curved beams 
cannot generalized quantitative terms. The response will vary according 
the general location the shear center, the relative magnitude the shear 


TABLE 3.—Shear Center Effect Dynamic Response Curved Beam Unsymmetric 
Cross Section 


Conditions 
(1) 
Variation (a) —4.975 in. 1.634 in. 33.5% stiffer 
Variation (b) 4.975 in. 1.634 in. 20.1% more flexible 


Note: in. 2.54 cm. 


center coordinates, and the shape and configuration the beam itself. was 
also found that for beams with small subtended angles, the shear center effect 
was significant completely alter the mode shape the dynamic response. 
The only conclusion possible that the shear center effects the dynamic 
response horizontally curved beams are significant, and should taken 
into consideration. 

Parametric Studies.—One the primary concerns this study was develop 
chart series charts which would readily applicable the determination 
the fundamental natural frequencies horizontally curved beams. order 
develop chart that would universal application, was necessary 
determine set parameters which would control the dynamic response 
the beam and define unique, correct solution the problem. 

result this study, the problem was defined according the following 
set parameters. The first parameter the ratio the torsional rigidity 
the beam the flexural rigidity the beam. This parameter was called the 
rigidity ratio and defined 


GK, 


a,= 


EM2 
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which was defined over range from 0.002 0.400. rigidity ratio 
0.002 roughly equivalent that isolated wide-flanged beam, while 
ratio 0.400 may considered (28) simulate the behavior curved 
beam built-up deck section, i.e., beam that supported and stiffened 
various means the web and flanges. should noted that the warping 
rigidity, has little effect the natural frequencies except for curved 
beams extremely short spans. 

second parameter slenderness ratio, defined L/r,, which 
length the curved beam, and polar radius gyration. The third parameter 
used was the subtended angle the curved beam. The problem was analyzed 
for range subtended angles from 90°. Finally, the problem was solved 
with variety end conditions. The three combinations end conditions 


10 20 30 40 50 60 70 80 90 
SUBTENDED ANGLE (DEGREES) 


FIG. 6.—Natural Frequencies Curved Beams Fixed-Fixed End Conditions, 
160 


studied were: (1) Both ends pinned; (2) both ends fixed; and (3) one end pinned-one 
end fixed supports. 

Having defined the foregoing controlling parameters the dynamic response 
the curved beam, numerous problems using different combinations the 
parameters were defined and solved the computer program developed for 
this study (9). The resulting solutions (eigenvalues) are generalized the charts 
the eigenvalue ratio /w,,, which fundamental frequency the 
curved beam under consideration, and w,, fundamental frequency straight 
beam having the same length, cross section, material properties, and end 
conditions the curved beam. 

selection the charts which resulted from this analysis shown Figs. 
through more complete range charts has been developed and presented 
Fehrenbach (9). 

Several general patterns dynamic behavior can readily observed the 
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charts. For very small subtended angles, the eigenvalue ratio approaches unity; 
that is, the curved beam behaves straight beam. the subtended angle 


SUBTENDED ANGLE (DEGREES) 


FIG. 7.—Natural Frequencies Curved Beams Pinned-Fixed End Conditions, 
160 


SUBTENDED ANGLE (DEGREES) 


FIG. 8.—Natural Frequencies Curved Beams Pinned-Pinned End Conditions, 
160 


increases, the curved beam becomes increasingly more flexible, and the eigenvalue 
ratio decreases accordingly. 
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CURVED GIRDERS 


Higher ranges denote stiffer torsional rigidities. torsional rigidities 
increase, the effect the beam curvature effectively diminished, and the 
curved beam behavior more closely approaches the equivalent straight beam 
behavior; accordingly, the eigenvalue ratio increases. Except for the low ranges 
torsional rigidity (a, 0.002) and slenderness ratio (L/r, 40), the effect 
the slenderness ratio the natural frequencies curved beams minimal 
(less than 10%). 

Finally, the end conditions have noticeable effect the curved beam 
behavior. Fixed-fixed end conditions contribute greatly the overall stiffness 
the beam, resulting higher eigenvalue ratios. can then expected, 
pinned-fixed end conditions are somewhat more flexible, with pinned-pinred 
ends being the most flexible combination end conditions. 


Summary 


Stiffness relationships have been formulated means variational procedures 
for the vibration problem spatially curved beams thin-walled open cross 
section. The resulting stiffness matrices are exact due the utilization shape 
functions which are derived from the solutions the homogeneous differential 
equations governing the static behavior. The mass matrices are general, that 
they include terms reflecting added inertial forces for singly symmetric 
nonsymmetric sections. course, for doubly symmetric sections, the mass 
matrices are identical the so-called consistent mass matrices. 

These element relationships were programmed for use digital computers, 
and solutions were obtained using the program that was developed. The results 
were then compared with existing solutions. series charts reflecting 
wide range applications was developed. Since the behavior curved beams 
coupled bending and torsion for the major axis, and bending and axial 
deformation for the minor axis, the structural response quite significantly 
affected the ratio flexural torsional rigidities. the charts, the parameters 
used are the ratio flexural torsional rigidities, the slenderness ratio with 
respect the polar radius gyration, the subtended angle the beam, and 
the type end conditions. The rigidity ratio ranges from 0.002, which represents 
most isolated wide-flange sections, 0.400, which reflects the equivalent rigidity 
ratio for curved girders deck system. has been noted that the actual 
torsional rigidity girder deck system much higher (tens and hundreds 
times) than that individual members (28). The slenderness ratio varies 
from 240, and the subtended angle from 90°. End conditions which 
were considered are pinned-pinned, pinned-fixed, and fixed-fixed. 

The charts can utilized effectively determining natural frequencies 
curved beams, demonstrated several examples. Good correlation was 
observed with existing solutions obtained Culver (6) and Shore and Chaudhuri 
(20). 

Examination the computed eigenvectors reveals that for doubly symmetric 
members with both ends pinned, the vibration mode shapes for translation and 
rotation are half sine waves, used most the classical solution techniques. 
However, for members with other than pinned-pinned end conditions, and shear 
center not coincident with the centroid, the mode shapes sometimes cannot 
defined simple trigonometric term. 
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addition, the element relationships developed herein have been utilized 
exclusively examples horizontally curved beams. However, the relationships 
deveioped herein can used describe the behavior curved beams 
any spatial configuration; for in-plane longitudinal vibrations. 
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cross-sectional area; 


coordinates shear center measured from the centroid; 
modulus; 
principal moments inertia; 
warping constant; 


stiffness matrix; 

St. Venant’s torsional constant; 
beam length; 

element length; 


T 


M,,M, moments about the and axes, respectively; 
mass matrix; 
[N] shape functions; 
axial force; 
{Q} generalized forcing function; 
amplitude beam vibration; 
polar radius gyration; 
T,,T, St. Venant’s and warping torsion; 


rigidity ratio; 
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kinematic degree freedom; 

displacement field; 

assumed displacement modes for dynamic behavior; 
vector nodal point displacements; 

variational operator; 

displacements the and directions, respectively; 
energy functional; 

mass density; 

kinematic degree freedom; 

angle rotation about the axis; 

angular frequency; 

curved beam natural frequency; and 

straight beam natural frequency. 
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INTRODUCTION 


There are many occasions when engineer has make judgment based 
his interpretation large amount rather imprecise information. One 
example the determination the cause causes structural failure 
and the synthesis general lessons can learn for the future. The final 
trigger incident which precipitates particular collapse may not the most 
interesting aspect it; very often there are complex underlying human interac- 
tions. Human error, one form another, has been shown dominant 
reason for many our failures (4,6,10). Indeed this conclusion should not 
surprise because the final analysis all error human error. people 
who have decide what do; people who have decide how should 
done; and people who have it. Engineers the past, flushed 
with the successes their applications the physical sciences engineering 
problems,have rather neglected their reliance human infallibility. engineering 
only the product, the hardware, physical system; the system which designs 
it, produces it, and uses it, human and therefore complex and vulnerable. 
There are many ways which the possibility human error the future, 
may reduced. For example, more attention might given the topic 
the education and training engineers; Melchers has pointed the need for 
better control the work (7). this paper, however, proposed concentrate 
one aspect the problem. How can conclusions distilled from past 
failures (and successes) when much the information available rather 
imprecise? The work described beginning; stems from the belief that 


there need for logical way examining information about failures, storing 
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them computer data structures and using these data basis for comparison 
with future proposals. The data structures are set using fuzzy logic. 
this paper attention will confined some the basic ideas fuzzy 
logic, with little mathematical, treatment, and illustration its application. 


INFORMATION 


Information can defined set propositions. The propositions may 
logically related and varying content and quality dependability. There 
not sufficient space here for full consideration these concepts (see Ref. 
but necessary outline briefly the difference between the dependability 
proposition and its truth even the probability its truth. proposition 
true then obviously dependable but the reverse not necessarily so. 
For example, modern relativistic physics has shown that Newton’s Laws are 
false, but nevertheless the predictions Newtonian mechanics embodied 
structural analytical techniques are highly dependable. modern philosophical 
view that all theories are models our understanding the world which 
are constantly attempting refute. matter how many times observations 
are made which agree with theory, not proved true; all that can 
said that highly tested, confirmed and dependable theory. the 
other hand, only one observation which disagrees with theory proves that 
theory false. This logical asymmetry between verification and falsification 
has been pointed out Popper (8). The argument stems from Hume who 
1739 first raised questions about the empirical explanation science 
process induction. showed that induction has assume that nature 
regular. Just because past futures have been like past pasts, does not 
logically follow that future futures will like future pasts. 

proposition precise, has high content and unlikely true, 
other words the probability that true low. the other hand, 
proposition vague, has low content and likely true, the 
probability that true high. For example, line measured and 
stated length 1.26359736 m’’ than highly likely that the 
last few decimal places are wrong. stated ‘‘the length between 
and then the proposition almost certainly true. This idea that 
there inverse relationship between the content proposition and its 
truth also comes from Popper (8). Engineers and scientists are not interested 
primarily theories and propositions which have high probability being 
true, but theories and propositions which have high content (and therefore 
low probability) but which are highly dependable. clear therefore that any 


measure the dependability proposition not the probability that 
true. 


Fuzzy 


Zadeh, who first published the idea fuzzy sets 1965, maintained that 
the way humans are able summarize masses information and then extract 
important items relevant problem, because think approximately. 
think terms classes sets objects where the transition from membership 
nonmembership not abrupt but gradual. example the fuzzy term 
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safe could defined fuzzy set the possible values probability 
failure, Fig. and term expensive could defined utility space, 
Fig. 1(b). Zadeh later suggested generalization ordinary two valued logic 
where the binary concepts true and false underlying the whole logic and 
mathematics, are generalized fuzzy sets the interval [0, 1]. 

The work described here based upon Baldwin’s considerable development 
this idea using the idea fuzzy truth restriction (2,3,5). Very little 
the underlying mathematics fuzzy logic will presented, instead attempt 
will made give some appreciation the ideas. 

Using the definitions made earlier may wish make the following simple 
proposition—Proposition: structure safe) then (it expensive) very 
true. This observation about the world. Suppose now particular observation 


1.0 

(a) (b) 

06 safe os expensive 

0.4 


FIG. 1.—Fuzzy Sets 


made quite safe what can deduce? The argument 
is, fact, logical modus ponens deduction the form: 


Proposition: structure safe implies structure expensive very true. 

Observation: Structure quite safe true 

Conclusion: The structure expensive recast the structure 
true. 


Expressed symbolically: ADB 


mathematics (5) allows the calculation and and are known. 
implication relation, such that Lukasiewicz, truth functionally 
modified The term the fuzzy membership function. The operation 

Another form logical deduction, which will used later, that modus 
tollens 
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Although fuzzy logic was developed generalization ordinary two valued 
logic, Baldwin recognizes that the interpretation the philosophical meaning 
the fuzzy truth restriction will depend the particular application. 
possible interpret measure plausibility, possibility, importance, 
dependability. this paper will regarded restriction upon dependability 
although the fuzzy logic notation will retained and will referred 
truth. Baldwin has defined standard fuzzy truth restrictions, some which 
will used the later example (Fig. 2). 


undecided 
unrestricted 


absolutely true Tabs 
very true 
true 

fairly true Trt 


€ 
® 
< 


FIG. 2.—Some Truth Restrictions, 


Before applying these ideas structural safety there one important 
characteristic the modus tollens deduction which must described and which 
relates Popper’s argument about the asymmetry between verification and 
falsification. may best understood simplifying the argument down 


that ordinary logic where each the propositions either true false. 
Consider the deduction 


Newton’s Laws (A) deflection beam (B) 
the beam deflects (B) 


Newton’s Laws are? 


The truth table ordinary logic Table 

Table are interest. However, this means that true false, 
that the observation does not verify nor falsify it, leaves undecided. 
Consequently fuzzy logic, the best (in the sense verification) that can 
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done the modus tollens deduction when and are fuzzy statements 
with fuzzy truth restrictions leave unrestricted. The best 
that can done falsify show that absolutely false. follows 


therefore that use modus tollens, the propositions should set that 
they may falsified. 


design perfectly safe) (S) then (it has low probability failures 
This means that necessary condition and sufficient condition 
NP. other words design safe then must have low probability 
failure but just because has low probability failure some defined 
limit state does not mean perfectly safe because other factors are also 
relevant safety. For example, design perfectly safe) (S) then (it has 
low probability failure due random hazards such fire, floods etc.) 


TABLE 1.—Truth Table Ordinary Logic for Implication 


Line number 


(1) 


(RH) and (it has low probability failure due human error) (HE) absolutely 
true. These may written: 


Note that includes the term perfectly safe which impossible satisfy. 
obtaining truth restrictions upon NP, RH, possible calculate 
the extent which false. other words will falsified truth restrictions 
upon NP, RH, and which represent the extent which believed that 
they falsify 

Now clearly NP, RH, are too general allow truth restrictions 


set them. clearly possible write down further necessary conditions 
for them follows: 


instruction and warning given user) WI] [(The structure not sensitive 
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(2) (4) 
true true true 
true false false 
false true true 
false false true 
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(No well-known mode failure missed) DEM 

(No calculation errors) DEE 

(Design consultants have adequate experience) DEC 

(Personnel available for design work are suitably experienced) DED 

(Personnel available for site supervision are suitably experienced) DES 

(Specifications are perfectly adequate) DEP 

[Quantity and quality research and development information 
sufficient] DER 

(There mode structural behavior inadequately understood 
existing technology) DET 

DEE (Checking procedures are adequate) DEEC 

DEE (Conditions work and employment are good) DEEW 

(Construction methods used are well-tried and tested) CEM 


(Construction company has adequate experience) CEC 
(Personnel available for falsework design are adequately experienced) 
CED 
[(Personnel available for site work are adequately experienced) CES] 
[(they are able appreciate technical problems associated with the design) 


CEP] 

CEE [(Record company good; mistakes infrequent past) CEEC] 
[(personnel working good conditions) CEEG] [(good happy relationship 
between personnel company) CEEH] (structure not sensitive erection 
method) CEEM 

(Perfectly normal contract procedure) CLC 

(No undue political pressures) CLP 

(No undue industrial pressures) CLI 

(No undue financial pressures) CLF 

(All random hazards have been identified) RHI 

[(The experience for this site suggests risk random hazards low) 
RHE (there change environment suggest this assessment wrong) 
RHC] [(perfectly adequate established data have been collected) RHD 
(the probability failure based these data low) RHP] 

{The notional probability failure [assuming perfect calculation 
model (LSM)] low} (NPF) 

(The probability new presently unknown effect the material 
causing failure low) (MMM) 

(The probability new presently unknown effect the structural 
form causing failure low) (SSS) 

(The probability new presently unknown loading effect causing 
failure low) (LLL) 

(The possibility that the calculation model perfect high) (LSMP) (the 
calculated notional probability failure perfectly dependable) 

(NPFR) 
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(The calculation model for the limit state perfect) LSM (the model used 
for the loads perfect) (LSML) 
LSM (The model used for the resistance perfect) (LSMR) 
LSML (The system used for the loads perfect) (LSMLS) 
LSML (The parameter statistics for the loads are perfect) (LSMLP) 


LSMLS (The model has been used before with problems) (F1) 


LSMLS (There change between this useage and previous ones) (F2) 


LSMLS (The test data available are perfectly satisfactory) (F3) 
LSMLP (The statistics have been used before with problems) (F4) 
abs 


LSMLP (There underlying change the useage) (F5) 


LSMLP (The data available perfect) (F6) 
LSMR (The system used for the resistance calculation perfect) (LSMRS) 
LSMR (The parameter statistics for the resistance calculation are perfect) 
(LSMRP) 
LSMRS (The system model has been used before with problems) (F7) 
LSMRS (There underlying change usage) (F8) 


LSMRS (Perfect prototype test data available) (F9) 
LSMRP (The parameter statistics for the resistance calculation have been 
used before with problems) (F10) 


LSMRP (There underlying change usage) (F11) 


LSMRP (Perfectly relevant measurements from sample surveys are available) 
(F12) 


These logical relations are summarized Fig. the figure the square 
boxes represent propositions and the circles are logical operators. Each proposi- 
tion has number (i) which its truth value Each logical operator may 
proposition may appear more than once the hierarchy, because its truth value 
stage the argument derives from the restriction the truth the 
proposition only one few the necessary conditions upon it. one 
proceeds the hierarchy from the input truth restrictions the bottom, following 
the arrows then more and more these restrictions are taken into account. 
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Finally, the restriction from the whole subhierarchy 
The intersection these three restrictions gives the final restriction 


The operations involving propositions with truth restrictions 41-46 are 
not described this paper (5). 


This hierarchy logical relations has been used obtain fuzzy safety measures 
for two famous failures, the Tay Bridge (9) and the Tacoma Narrows Bridge 


TABLE 2.—Assessmenis for Examples 


Truth Restriction 


absolutely false very false false 
very false absolutely fairly false 

absolutely false absolutely false unrestricted very false 

absolutely false false unrestricted fairly false 

very false unrestricted fairly false 

absolutely false unrestricted fairly false 

false unrestricted unrestricted 

absolutely false unrestricted very false 

absolutely false absolutely false unrestricted false 

very false unrestricted 

absolutely false 

absolutely false unrestricted fairly false 
unrestricted unrestricted 

unrestricted unrestricted very false 

absolutely false unrestricted false 

very false absolutely false fairly false false 
unrestricted unrestricted 

unrestricted absolutely false unrestricted unrestricted 

unrestricted unrestricted unrestricted fairly false 

unrestricted fairly false unrestricted very false 

unrestricted unrestricted 

false unrestricted fairly false 
unrestricted absolutely false 


unrestricted unrestricted absolutely false 
unrestricted 


unrestricted 
unrestricted 


(1). The safety measure fuzzy truth restriction which falsifies Subjective 
assessments for the input propositions the bottom the hierarchy are 
shown Table computer program available process these inputs 
and produce fuzzy truth restrictions every proposition every stage 
the hierarchy. 


The resulting truth restrictions for the Tay Bridge assessments are shown 


EM2 
sition Tacoma Tay sition Tacoma Tay 
(1) (2) (3) (4) (5) (6) 
F9 
F10 
MMM 
SSS 
LLL 
RHD 
RHP 
RHE 
RCH 
RHI 
CLP 
CLI fairly false 
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Fig. The numbers the truth restrictions that figure must related 
the propositions the boxes Fig. The critical truth restriction which 
yields the most false restriction upon which derives from which 
turn derives from the restriction absolutely false given the proposition 
LLL (that the probability new presently unknown loading effect causing 
failure low). The truth restrictions deriving from and are and 
Tespectively and have similar values. The truth restrictions from which 
derive are shown Fig. 4(b). These show that the design errors the 
severity the the construction errors and the misuse 
the structure are rather similar. The restriction upon derives 


FIG. 4.—Results for Tay Bridge 


FIG. 5.—Results for Tacoma Narrows Bridge 


directly from the restriction absolutely false RHD which refers the lack 
data about the severity wind storms the Tay Estuary. 

The truth restrictions for the Tacoma Narrows Bridge assessments are shown 
Fig. This time and for and are identical and for 
undecided. Here derives directly from the assessment absolutely false 
SSS which refers the new mode oscillation the bridge due its unusual 
structural form. The restriction upon comes from with little 
reason falsify from construction error, use climate (CE,SU,CL). The 
restriction derives directly from the input absolutely false DET 
the proposition that there was mode structural behavior inadequately 
understood the existing technology the time. 


1 
(a) 
Xt ) 
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The inputs Table seem therefore produce truth restrictions 
NP, HE, and which agree with what might expected intuitively from 
the mass data about the two accidents. 


Future Use 


There are two major problems which must dealt with this analysis. 
Firstly could argued that the truth restriction upon for structure which 
know have failed should absolutely false. Secondly, the answer derives 
from assessments made with the benefit hindsight and therefore little 
use for comparison with data regarding the future. order correct both 
these difficulties necessary for research such that Walker and 
Sibly (10), which attempts understand the thoughts of, and the pressures 
on, the designers these bridges. Just what assessments might they have made? 
The reason that the fuzzy truth restrictions are not absolutely false 
that not possible have totally clear idea the cause failure and 
every truth restriction the critical path through the logical hierarchy 
not absolutely true false. With increasing clarity certainty regarding the 
reasons for the actual failure can the uncertainty reduced. 
The restrictions Figs. and contain absolutely false. The truth modifiers 
contained within the hierarchy Fig. have empirically adjusted 
reflect the available information and the result stored computer data file. 
anticipated that this process could performed all failures from the 
published information where incomplete information available then the truth 


restrictions are simply left unrestricted. Work proceed methods 
combining these files sort learning procedure which would provide 
accumulated experience within the computer against which all new assessments 
may compared. 
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ULTIMATE BEHAVIOR I-BEAMS 
Yuhshi ASCE and Susumu 


INTRODUCTION 


The purpose this study derive fundamental equations single 
curved I-beam subjected the action bending and torsional moments and 
also investigate the behavior curved flexural members under large torsional 
deflection. The derived equilibrium equations are solved the transfer matrix 
method for both elastic and inelastic ranges. 

Recently, the Task Committee Curved Girders (3) presented design guide 
for curved steel I-girder bridge and the Task Committee Curved Box Girder 
(4) summarized the researches conducted investigate the static and dynamic 
behavior curved steel box girder bridges. 

The two committees have identified one problem among others—relatively 
little data are available for use developing ultimate load design criteria for 
curved girders and the committees suggested that further research ultimate 
load behavior needed for analytical procedures. Keeping the foregoing points 
mind, the writers have proceeded with the present study. 

Experimental studies are also made determine the effect curvature 
the ultimate load behavior. Six simply supported single curved I-beams are 
tested under concentrated load. Curvatures the I-beams vary from the 
straight beam with specified initial crookedness, f,/L 1/1000, 
the curved beam L/8R 1/100, which the curved length 
the beam; and the radius horizontal curvature. 


Assumptions.—The following assumptions are made the analysis single 
horizontally curved I-beam shown Fig. (1) Length and radius curvature 
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are very large compared with the cross-sectional dimensions; (2) cross-sec- 
tional distortions occur during deflection; (3) shear strain along the middle surface 
the thin walls due shear equilibrium with the change normal stress 
small and may neglected; (4) shear strain the plane normal the 
middle surface the thin walls small and may neglected; and (5) transverse 
displacements are much larger than the longitudinal displacements. 
Strain-Displacement Relations.—As shown Fig. the x-axis and y-axis 
coincide with the respective principal centroidal axes the cross-section and 
the z-axis tangent the curved axis the member. The origin the cylindrical 
coordinates (p, and taken the center curvature. 


FIG. 1.—Analytical Model Curved Member 


From the finite displacement theory (11) the following strain-displacement 
relations may obtained considering the fifth assumption 


which and denote the displacements the and directions, 
respectively. The second assumption leads the following equations: 
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where and the displacements the centroid from the original position; 
and the angle rotation about the 
Using Eq. and the third and fourth assumptions which imply y,, 


the middle surface, the longitudinal displacement, may derived 
shown: 


which the terms with prime denote the partial derivatives with respect 
and constant (i.e., with respect and the longitudinal 
displacement the centroid. The integral terms Eq. may transformed 
using the curvilinear coordinate axes which taken along the middle 
surface the thin walls; and normal the s-axis. Then, Eq. may 
rewritten the following convenient form: 


which denotes unit warping with respect the centroid and may simply 
defined the following equation for the approximation R/p* 


The terms with asterisk, p*, x*, and show the values respectively, 
the middle surface the thin walls. 

The general equation shear strains any surface parallel the middle 
surface the thin wall may given 
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With the aid Eqs. and Eq. may formulated shown: 


which the higher-order terms are small and are neglected. 


Substituting Eqs. and into Eq. and considering the first assumption, 
R/p the following equation may deduced 


Equilibrium Equations and Boundary Conditions.—The curved member this 
analysis equilibrium under the action distributed loads along the length 
and external forces both ends. Using the variational method (8-10), internal 
and external works may formulated 


L 


L 
4 


L 


end forces both Substituting and into Eq. and 
requiring that the equation satisfied for arbitrary virtual displacements leads 
the equilibrium equations and the boundary conditions, follows: 


” 


(6c) 
L 
Wo 
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A A A 


A A 


A 


The force-displacement relations for doubly symmetric cross sections may 
obtained substituting Eq. into Eq. shown 


sin 


A 
A 
A 


€,0, wo R 
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equal the normalized unit warping for doubly symmetric cross sections. 
Eq. lle represents St. Venant’s torsion which naturally appears because 
the variation shear strain through the thickness the plate. 

The procedure followed the derivation Eqs. 9-11 similar that for 
the straight members described Refs. and except the necessary changes 
required for the curved members. 

Transfer Matrix Method.—From the boundary conditions (Eqs. 9e-9h) the 
displacements and the forces the state vectors may derived shown: 

Wo , Yo M M’ K 


which the terms with the denote the bending moments with respect 
global coordinate system. 
The equilibrium equations from Eqs. 9-12 may represented matrix 


The field transfer matrix relating the state vectors the left side the 
right side the divided elements may derived integrating Eq. using 
the Runge-Kutta Method. The linear analysis curved beams the transfer 
matrix has been presented Ref. computer program was developed here 
provide the numerical results utilizing the successive approximation method 
for the solution the nonlinear matrix equation (Eq. 13). 

Stiffness Partially Yielded Cross Section.—Using the force-displacement 
relations for the doubly symmetric cross section, Eq. may transformed 
into: 


=x 
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which denotes the polar radius gyration area. the foregoing equation, 
the last term small and may neglected. that case, may noted 
that changes linearly with respect and w,. However, does not 
vary linearly with respect displacement quantities, may seen from 
Egs. and 14. 

the inelastic range, penetration yielding reduces the stiffnesses the 
member. Eq. applicable the partially yielded cross section where the 
stiffnesses change from elastic inelastic ones. For the ideal elastic-plastic 
stress-strain relationship, the inelastic stiffnesses the cross section may 
derived using the secant modulus elasticity shown 


A A 
The inelastic stiffnesses are, thus, determined trial and error procedure using 


and Using the total strain theory, St. Venant’s torsional stiffness 
the partially yielded cross section may derived follows 


Ap 


which and correspond the elastic and plastic parts, respectively. 


Six welded curved I-beams are tested investigate the load-deflection relations 
and the ultimate load behavior under concentrated load the span center. 

Curved I-Beams and Test Rig.—The cross-sectional dimensions and the radius 
curvature test beams are summarized Table Beams AR-3 and BR-3 
are designated the laterally unsupported straight beams with the initial 
crookedness L/1,000. The curved flange plates are cut out from 8-mm thick 


TABLE 1.—Dimensions and Curvatures Curved 


Length Span along the Curved Web Center Line 
Nominal Value and Measured Value 
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plate, and the flat web plate bent elastically and welded into the position 
give specified curvature. 

these experiments, beams are restrained their supports described 
Ref. providing restraint against twist but not against warping. Load 
applied the center the compression flange through the sole plate 
tension jack through the loading frame shown Fig. 
The jack which attached the load cell then connected the gravity-load- 
simulator hold the applied load vertical during test. 


FIG. 2.—Test Setup 


Initial Imperfections.—The longitudinal residual stress distributions obtained 
the sectioning method for AR-2, and AR-3 beams are presented 
Fig. The solid line curve through the experimental results considered 
the representative distribution curve and may defined the following 
equations: 
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meganewtons per square meter; 


meganewtons per square meter 


The foregoing equations satisfy the self equilibrium conditions resultant axial 
force and bending moments about both axes shown: 


A A 


Measurement out-of-straightness were made seven equally spaced points 
for each fabricated curved specimen. Fig. shows v,, and curves for 


AR-I (17100) 
AR-2 (1/200) 


o—-—o 


4 


FIG. 3.—Measured Residual Stress Distributions 


the BR-series. The curved solid thin lines indicate the prescribed curves 
for the test beams. found that there less than 0.7 difference between 
the prescribed and the obtained curves. 

Ultimate Load Behavior.—Fig. shows the relations between the load and 
the midspan displacements and for BR-1. The test results are represented 
the thick lines joining the test points and the thin curved lines represent 
the numerical results described the next section. load increased, 


300 
200 
200 
-200 
Va 
200 
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the beam stressed elastically until parts the beam are strained into the 
inelastic state. The ultimate capacity reached when the progressive yielding 
reduces the stiffness the beam zero. Good agreement may seen between 
the test and numerical results for the beam. 

Figs. and show the relationship between midspan load and angle rotation 
midspan for the AR-series and BR-series, respectively. The test results are 
represented the solid lines joining the points marked with small circles, 


R Uo 


1 mm 
Radian 


FIG. 4.—Shapes Geometrical (BR-Series) 


and the broken lines depict the behavior predicted the present study. 
most the cases, the experimental curves and the theoretical curves are very 
close. the aforementioned figures, seen that the load-deformation relations 
for AR-3 and BR-3 correspond the lateral buckling behavior beams with 
initial crookedness The inset Figs. and lists the experimental 
and theoretical ultimate loads the test beams. 


Elastic Analysis.—An example the numerical results obtained under the 
loading equal end moments shown Fig. simply supported beam 


BR-3 
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(2.1 10° and 0.3 used for numerical computation. Parameter 
L/8R which similar rise/span ratio circular arches introduced here 
for the horizontally curved configuration. 

The curved length the beam Fig. corresponds 1.0 which 
and the modified beam slenderness (7); and the elastic 
critical moment laterally unsupported straight beam with the span length 

Fig. the curved thick lines show the elastic large-deformation behavior 
predicted the present study, and the thin lines show the simplified results 


obtained considering only the first term Taylor’s expansion whereby 
Egs. 9-11 are linearized. 


ton 


FIG. 5.—Load Versus Deflection Relationship (BR-1) 


Fig. the critical elastic moment (10) curved beam and 
defined as, 


cr/ App R? L? T L? 


Eq. determined the moment when the displacement (the curved 
thin lines Fig. approaches infinity. interesting note that the u/b-value 
increases gradually with M,/M, until the u/b versus M,/M, curve reverses 
the direction the elastic large deflection analysis while the 
and increase with the value 


Ton 
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—o—o— Test results 
Numerical results 


0.02 003 0.05 007 Rad. 


FIG. 6.—Load Versus Torsional Angle Relationship (AR-Series) 


Test results 
Numerical results 


FIG. Versus Torsional Angle Relationship (BR-Series) 
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Inelastic Analysis.—Fig. shows the results the inelastic analysis the 
same beam Fig. The curved thick lines are for the inelastic analysis, 
while the thin lines are for the elastic analysis. The figures the left show 


FIG. 8.—Elastic Analysis Curved Beam 1/100) 


FIG. Analysis Curved Beam 1/100) 


the inelastic behavior and the resulting ultimate moment without residual stresses 
and the figures the right show the same with specified welded-type residual 
stress distribution shown the inset Fig. seen from Fig. that 


(Mcrlapp 
0.50 
025 050 075 1.00 
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Case (a) 
83Mp 
(test) (b) 


Case (c) 


Theoretical Result 
Case|Mork| C | 


Ultimate Strength 
L 284.8cm 


FIG. 11.—Approximate Ultimate Strength 
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the early yielding due the residual stresses reduces the ultimate moment 
from 0.53 M,-0.49 M,. 

The numerical results obtained for the AR-series and BR-series are shown 
together with the test results the inset tables Figs. 5-7. The modified 
residual stress distributions expressed Eq. and the prescribed curves 
are used the inelastic analysis and for the ultimate strength calculation. 

Ultimate Strength for Certain Loading Conditions.—The ultimate strength for 
the loading conditions (1) Equal end moments; (2) uniformly distributed load; 
and (3) concentrated load the midspan determined separately for the 
three cases. Fig. shows the relation between the obtained ultimate strength 
and corresponding L/8R 1/200 and 1/50. When 
load applied the centroidal axis the beam, the may obtained 
using the following values (2): 


The inset Fig. shows the residual stress pattern used this analysis. 

Fig. 10, the calculated ultimate moments marked triangle and circles 
for the three different loadings are nearly the same for the same and 
values. The thick solid lines represent the approximate ultimate moment calculated 
the formula covered the following section. 

Approximate Ultimate Strength.—Approximate ultimate strength may simply 
determined the point intersection between the second-order elastic 
curve and the second-order rigid plastic curve shown Fig. 11. The dimensions 
the test specimen are given Fig. 11. The chain curves marked 
Matrix represent second-order elastic results obtained from the numeri- 
cal solution Eqs. 9-11. approximate second-order elastic curve, leading 
relationship for the ultimate strength, can developed follows. 

The horizontally curved beam could defined cartesian coordinate system 
with the z-axis connecting the ends the beam. Then the initially circular 
horizontal curve approximate 


and also assumed sinusoidal the form 


The displacement components and are measured from the initial curved 
configuration, that the total displacement measured from the z-axis 


and the total rotation 


The differential equations with initial deflection and under equal end 
moment are given (12): 
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with the following solutions for and 


2 
ro Pr 


Now, defined the total lateral displacement the compression flange 
midspan. Therefore, due the horizontal curvature, follows that 


+ 


“,= 
tf 
2M, 


seen from the figure, the curves resulting from Eq. closely coincide 
with the second-order elastic curves obtained for the curved beam for L/8R 
1/100 well for 1/1,000. 

The dropping curve represents the second-order rigid plastic curve (6) which 
governed the lateral direction plastic strength the compression flange 
bending beam. This relation can developed follows. 

Under the axial compressive force the flange, plastic moment due 
the lateral bending the compression flange may expressed 


M,b 


Considering M,, u,, Eq. may reduced the following Eq. 
which the equation the dropping curve Fig. 


382 EM2 
With simply supported boundary conditions, possible solution given 
(24a) 
M, M, 
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Approximate ultimate strength M,/M,) the point intersection 
between the second-order elastic curve and the second-order rigid plastic curve 


may determined from the following resulting equation obtained equating 
and 27. 


2M, \2Rb 


The approximate ultimate strength 0.439 and 0.792 shown Fig. 
are then obtained from Eq. for L/8R 1/100 and 1/1,000, respectively. 
The approximate ultimate strengths for L/8R 1/200 and 1/50 are shown 
the thick lines Fig. 10. From Fig. 10, seen that there good 
agreement between the theoretical ultimate strength results (marked 
and the approximate strength curves and may noted that the approximate 
curves are the conservative side compared the theoretical results. 
Two results for AR-2 (test) and BR-2 (test) and marked such are also 
shown for comparison with the calculated results. 


method analyzing the elastic and inelastic large-deflection behavior 
the horizontally curved I-beam presented for certain loading conditions using 
the transfer matrix method. Load versus and relations are investigated 
into the elastic finite displacement range shown Fig. and the results 
the inelastic analysis considering the residual stress distributions are shown 
Fig. From Fig. 10, seen that the ultimate strength for the three 
different loading conditions are almost the same. 

The equilibrium equations are then reduced the linear simultaneous differen- 
tial equations and the critical elastic moment thus obtained. Approximate 
ultimate strength determined the point intersection between the second- 
order elastic curve which asymptotic the bifurcation solution and the 
second-order rigid plastic curve. Good agreement found between the theoretical 
and the approximate results. 

The experimental studies for six beams also show good agreement with the 
theoretical ones respect the and curves and the ultimate strength. 
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The following symbols are used this paper: 


area cross section; 
modulus elasticity and secant modulus; 
shear modulus; 
moment inertia about x-axis and y-axis; 
warping moment inertia; 
St. Venant’s torsional constant; 
cross-sectional value; 
length curved I-beam; 
approximate critical moment curved beam; 
elastic lateral buckling moment; 
full plastic moment; 
moment about x-axis and y-axis defined local coordinate 
system; 
moment about x-axis and y-axis defined global coordinate 
system; 
warping torsional moment; 
axial force; 
radius curved beam; 
St. Venant torsional moment; 
curvilinear coordinate axes; 


(p, 
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coordinate axis connecting ends beam; 
shear strain; 

nondimensional moment; 

strain; 

modified beam slenderness; 

tangential and normal components distance between point 
cross section and centroid; 

cylinderical coordinate axes; 

normal stress; 

angle torsion; 

unit warping with respect centroid; and 
normalized unit warping. 
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ELASTIC STRESSES FRETTING FATIGUE 


INTRODUCTION 


Fretting fatigue the combined action fretting and fatigue. Fretting corrosion, 
which the result minute relative movements between two bodies that are 
held contact normal force, was first examined detail Eden, Rose, 
and Cunningham 1911 (2). Extensive experiments have been reported since, 
and various theories have been advanced explain the fretting mechanism. 
For brief survey the literature the subject, the reader may wish 
refer Refs. and 10. 

some the recent fretting fatigue experiments conducted the third 
writer (who recorded his observations unpublished report), cleavage cracking 
was observed immediately below the contact area. possible explanation 
this observation has been suggested that high triaxial strain developed 
below the contact surface which, the presence some localized discontinuities 
(such micropores included particles), could result localized cleavage 
fracture. 

light the foregoing observations, appears desirable take close 
analytical look the distribution stresses the vicinity the contact area 
typical fretting-fatigue problem. The specific physical problem considered 
herein for this purpose (see Fig. consists large rectangular plate 
unit thickness, and two smaller fretting pads, also unit thickness, which 
exert pair equal and opposite transverse forces, shown the figure. 
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The plate fixed along one end and subjected cyclic axial force, 
acting the other end. Under the action the cyclic load, the plate elongates 
and shortens function time, and frictional forces develop between the 
plate and the fretting pads. The specific purpose the present study 
analyze the stresses produced the plate the combined action the cyclic 
axial load, and the normal and frictional forces that occur between the plate 
and the fretting pads. 


The analysis presented herein based the following assumptions: 


The material homogeneous, isotropic, and linearly elastic. 

The frequency the cyclic load much lower than the lowest natural 
frequency the plate, thus implying that the problem does not involve any 
substantial dynamic effects. 

The plate thin, that the problem can approximated 
problem. 

The length the contact area much smaller than the dimensions 


FRETTING PAD 


2b<<Ho 


Q=CYCLIC LOAD 


FIG. 1.—Physical System Analyzed 


particular, that analyzing the stresses the vicinity the contact area, 
the finite plate can approximated semi-infinite one, thus making possible 


take advantage some basic analytical that are already available 
the literature (8). 


Letting and denote the stress exerted the fretting 
pad (see Fig. one can write 


FRETTING FATIGUE 


which friction coefficient; and outstanding height the fretting 
pad. Note that Eq. based the assumption that the fretting pad behaves 
essentially column fixed against rotation both ends. 


MAXIMUM LOAD 


FIG. 2.—Loads Used Analysis 


FIG. 3.—Distribution Contact Loads: (a) Effect (b) Effect 


The actual distribution the contact forces not known. For analytical 
convenience assumed that the distributed normal loads, associated 
with and vary accordance with the ordinary theory beams, i.e. 
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which moment inertia; the coordinate system (x,y) shown 
Fig. and and q,, are considered positive when compressive. further 


assumed that separation occurs between the plate and the fretting pads, 
i.e. 


all contact points. Since slippage considered have taken place 
few cycles loading, appears reasonable express the friction force 


FIG. 4.—Classical Solutions: Normal Force; Tangential Force 


q (x) =ARBITRARY DISTRIBUTED LOAD 


FIG. 5.—Fixed and Moving Coordinate Systems 


distribution the following. [Analytical results based alternate 
assumption (namely, the assumption parabolic distribution that shear 
stresses the ordinary theory beams) can found Ref. 1.] 


except near the corners the fretting pads where from theoretical 
considerations. 


Two 


the special case semi-infinite plate subjected concentrated normal 
load, shown Fig. 4(a), the stress distribution given (8) 


| 
y 7 
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the case tangential force, shown Fig. 4(b), the corresponding 
expressions are: 


respectively. These two classical solutions, used conjunction with the principle 


superposition, form the basis all the analytical results presented this 
study. 


the two cartesian coordinate systems shown Fig. (x,y) fixed, 
and moving, subject the conditions 


which relative distance between the origins the two systems. 
Consider now distributed normal load, q(x), acting the edge the 

semi-infinite plate, shown Fig. The stresses produced can 

obtained from and substituting g(c)dc for and integrating. Thus 


EM2 391 
T r 
or, cartesian coordinates 
cos 


APRIL 1981 


Stresses Produced P.—The distributed contact load, associated with 


given Eq. Substituting Eq. into Eqs. and integrating from 

(o,)p 2 2 2 2 


xy? 


Stresses Produced M.—Similarly, the stresses produced can obtained 
simply substituting Eq. into Eq. and integrating. Thus 


which log denotes natural logarithm. 
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Stresses Produced V.—By analogy 12, the stresses produced 
distributed tangential load, s(x), can obtained from Eqs. and 
substituting s(c)dc for and integrating. Inasmuch s(x) associated with 
the present problem, one can write: 


Rewriting Eq. 


| 


separate the contributions and can expressed as: 
(o,) ve + vm 
which, after some lengthly mathematical manipulations 
2 2 2 2 + log 2 2 


tan 
n um 2 y 


2 


tan 


(17a) 

(185) 
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Summary Stresses.—The stresses produced the plate and 
decrease rapidly the distance from the contact area increases. follows 
that, analyzing the stresses the general vicinity particular fretting 
pad, the contribution the other fretting pad can neglected without introducing 
any substantial errors. The overall stresses arbitrary point interest 
can thus expressed as; 


«@,,) VM 


which 


represent the stresses produced the maximum tensile value the cyclic 
load, 


was observed various investigators (3,6,7), the friction coefficient, 
appears vary with the number cycles during fretting test. the 
present study, numerical solutions were obtained for the following specific cases: 
0.2, 0.4, 0.6, and 0.8. The other major problem parameter that was 
varied this study 


which (Q,/H,) maximum value axial cyclic stresses; and (P/b) 
twice the average normal contact load under the fretting pad. Numerical solutions 
were obtained for and 10. The fretting pad height, was assumed 
have fixed value, all cases that involved axial load, 
were considered moving), the ratio was treated variable problem 
parameter, and numerical solutions were obtained for several different values 
h/b the range 0.2 h/b 3.0. (Note that this connection cannot 
exceed since otherwise the condition separation, i.e., Eq. would 
not satisfied.) 


(21a) 
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For each set values the problem parameters, the state stress was 
determined numerically total number 861 points the general vicinity 
each fretting pad, namely, all points grid system defined Ax, 


(a) 


(c) Tmax 


FIG. 6.—Stress Contour Lines 0.4; h/b 0.2) 


and the maximum shear stress were then obtained each point treating 
the problem two-dimensional one, i.e., from the equations 
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Note that the third principal stress this plane stress problem. 
When o,,,, and Eqs. are different signs, Eq. does indeed 


(a) max 


(b) omin/ 


(c) Tmax 
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define the maximum shear stress three-dimensional sense. When and 
have the same sign, however, would have calculated one 
half the absolute value Whichever happens numerically 
larger. 


-0.55 


(c) Tmax 


effort visually enhance the significance the numerical results 
obtained this study, stress contour diagrams were constructed for 
and for each set values the problem parameters considered. Due 
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space limitations, however, only small number these graphs are presented 
this paper. 


(a) 


(b) min 


(c) Tmax 


Fig. corresponds relatively large value the axial cyclic stress 
5), whereas Figs. 7-9 are for relatively small value that stress 
1), with varying the range 0-0.4. The numerical coefficients shown 
these four figures are multiplied determine the actual 
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stresses. effort illustrate the distribution strains within the plate, 
the results corresponding one particular case 0.4), and obtained 
using value 0.34 for Poisson’s ratio are summarized Fig. 10. 


FIG. 10.—Stress Contour Lines 0.4; h/b 0.2) 


determine the actual strains. Finally, Figs. and are for the 
special case axial cyclic load) with h/b 0.2 and 
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0.2, 0.4, respectively. Since these two figures, the quantity 
can longer used common multiplying factor. Accordingly, 


the quantity used for this purpose. 


-0.44 


(a) oT max/ (P/b) 
-0.70 


(b) 


Tmax 


seen from examination Figs. that the stresses produced 
the contact loads and represent essentially localized phenomenon, 
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cases considered, the maximum tensile stress occurs the edge the plate 
immediately outside the contact area, whereas the location the maximum 
shear stress varies substantially from one case another. most cases, the 


0.30 =>+0.42 


(a) 


FIG. 12.—Stress Contour Lines 0.4; h/b 0.2) 


shear stress reaches local maxima several different points, some which 
lie the surface and some below it. Fig. e.g., there are total number 
five local maxima for the shear stress, the absolute maximum among them 
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occurring directly under the fretting pad depth about important 
recall this connection that Figs. 6-12 refers maximum shear 
stress the plane stress sense. When and have the same sign, 
three-dimensional sense would have obtained one half the absolute 
value the numerically larger principal stress. 

According the linearly elastic analysis presented herein, the question 
crack initiation and propagation would normally most critical the edge 
the plate immediately outside the contact area since: (1) The maximum values 
the principal stresses and strains occur that region; and (2) the surface 
abrasion associated with the fretting action would conceivably produce additional 
stress concentrations not included the present analysis. This observation 
consistent with experimental results reported various investigators (4,9). 

conceivable, the other hand, that cracks could also initiate some 
subsurface points result various contributing factors, such as: (1) Stress 
concentrations associated with some localized discontinuities; (2) deviations 
from the linearly elastic stress distribution presented herein due residual 
stresses other nonlinear behavior. 


Summary 


The distribution stresses typical two-dimensional fretting-fatigue problem 
has been analyzed assuming homogeneous, isotropic, and linearly elastic 
material. solutions have been obtained that are within 
the assumptions used. Although extensive numerical results have been obtained 
during the course the investigation, only limited results have been included 
this paper due mainly space limitations. Inasmuch the solutions presented 
are essentially exact and the closed type, the interested reader could generate 
numerical data for any specific set values the problem parameters with 
minimal effort. 

According the numerical results presented herein, small region the 
plate immediately outside the contact area would appear most critical 
from the viewpoint crack initiation. However, under certain conditions (e.g., 
the presence micropores other types discontinuities) cracks could 
also initiate some subsurface sites. 
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The following symbols are used this paper: 


one half contact area length; 
distance between fixed and moving coordinate systems (see Eq. 11); 
depth plate; 

height fretting pad; 

length plate; 

Vh/2 moment; 

normal force exerted fretting pad; 
maximum value cyclic load; 
distributed load due 
distributed load due 

distributed load due 

tangential force; 

frictional coefficient; 

Poisson’s ratio; and 


Q./H,. 


h = 
s(x) 
= 
v= 
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ANALYSIS UNBALANCED 
ANISOTROPIC SANDWICH PLATES 


INTRODUCTION 


Rigorous solutions for certain types sandwich plates with unbalanced 
laminated faces have not been realized due the existence anisotropic coupling 
stiffnesses. Typical examples such plates are those with unbalanced angle-ply 
faces. Neglecting the coupling action can result appreciable errors the 
deflection, stress, and frequency predictions. 

Recently, modified stiffness formulation was presented Ref. which 
allows the uncoupling the bending and the membrane actions sandwich 
plates. result, the analysis unbalanced anisotropic sandwich plates under 
static loading that were not amenable rigorous solutions has been realized. 

The objectives this paper are extend the modified stiffness method 
the dynamic analysis anisotropic sandwich plates for which rigorous solutions 
are not available, and demonstrate the effect bending-membrane coupling 
the static and dynamic structural responses sandwich plates with unbalanced 
angle-ply laminated faces. 

The modified formulation presented herein independent the boundary 
conditions the plate and can used conjunction with variety analysis 
techniques. order demonstrate the potential the formulation, solutions 
are generated herein for simply-supported sandwich plates using double Fourier 
series procedure. The modified formulation yields solutions that are complete 
agreement with those using the formulations presented Refs. and for 
sandwich plates with unbalanced cross-ply faces. The method, therefore, should 
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accurate for other sandwich plates, such those with unbalanced angle-ply 
faces. 


Generat 


The sandwich plate shown Fig. consists orthotropic core with 
shear moduli and G,, and two unbalanced laminated face plates subscribing 
the usual thin plate assumptions which transverse shear deformations are 
neglected. assumed that the core incompressible the transverse direction 
and carries longitudinal stresses, but that transverse shear strains are 
significant. The two unbalanced anisotropic laminated face plates are assumed 
have the same thickness, number plies, and material properties. 

The uncoupled strain energy the unbalanced anisotropic laminated face 
plates, U,, can expressed terms the modified stiffness matrices [R*], 


FIG. Plates and Displacements 


[D*], and the partial curvatures and (1): 


ax? axay 


which 


Also, (w,,,W,,) and are the partial deflections due bending and 
transverse shear deformations, the and the directions, respectively. 


Kinetic 


The total kinetic energy the sandwich plate considered the sum the 
kinetic energies the core, and the face plates, 


é 


DYNAMIC ANALYSIS 


Based the assumptions considered, that the faces the sandwich plates 
have the same thickness, number plies, and material properties, the centers 
the mass the sandwich plate the and directions are located 
the middle plane the core, which turn, coincides with the reference plane 
the sandwich plate. Also, the centers the mass each face the 
and the directions are located the reference plane the individual face 
plate, since the mass density ply assumed independent the 
orientation the ply. result, for the sandwich plates considered the 
work, the bending and membrane actions are not coupled the kinetic energy 
expressions the face plates and the core. 

For the case transverse vibration, the kinetic energy the two face plates 


that includes both the translatory and rotary inertia terms can expressed 
follows: 


and the mass density the face plate. 


The kinetic energy the core that includes the rotary inertia terms can 
expressed 


the preceding equation 


ze 


which the mass density the core. 

Note that both the modified strain energy expression (Eq. and the kinetic 
energy expression (Eq. are functions the partial deflections only. The 
governing differential equations motion generated from the application 
the variational principle can expressed the following: 
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the preceding equations, and the transverse shear forces the 
orthotropic core (1); and 


a’ (w,, 
a? (W,, + W 


= —a,(W,, + +a, 2 


and 


pointed out that the difficulty associated with the anisotropic coupling 
the governing differential equations (Eq. has been totally 
eliminated, since addition the uncoupling the energy expressions, all 
the anisotropic stiffnesses vanish Eq. (1). 


Note that the formulation independent the boundary conditions the 
plate, and that the governing differential equations are function the partial 
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deflection only. such, variety analysis procedures may applied. 
Herein, order demonstrate the potential, and assess the accuracy the 
formulation, Fourier series approach applied the prediction the dynamic 
response rectangular simply-supported sandwich plates. 


The boundary conditions for rectangular simply-supported plate 
are follows (1): 


for y=0,B 


The boundary conditions (Eq. 13) and the differential equations (Eq. are 
satisfied the following series expansion: 


(Wo, Wy Wax) (W W oymn ’ sin L sin B (14) 


m=i n=l 


and into the governing differential equations (Eq. yields the following 
system equations: 


(2c) Angle-ply faces 


FIG. 2.—Sandwich Plate and Ply Lay-Up Configuration 
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sxmn 


Solutions Eq. yield the frequencies w,,, the structure for each 


The analysis this section focused demonstrating the potential, and 
assessing the accuracy the formulation predicting the frequencies sandwich 


TABLE 1.—Properties Face Laminates 


Density, 
pounds- 
Face square seconds pounds per 
material per square inch 


(1) (2) (3) (4) (5) 


0.15 7.5 10° 


TABLE 2.—Properties Cores 


G,, G,, 

Density, pounds pounds 
pounds- per per 

Core square seconds| square square 
designation Core material per inch inch 


(1) (2) (3) (4) (5) 


Cellular cellulose 0.22 500 


plates with unbalanced anisotropic laminated faces such those with angle-ply 
faces for which rigorous solutions are unavailable, and illustrate the effect 
coupling the static and dynamic structural responses sandwich plates 
with angle-ply faces. 

pointed out that coupling between membrane and bending actions 
one the main characteristics unbalanced laminated constructions. For certain 
types sandwich plates with unbalanced laminated faces [see Fig. 2(b), 
and Fig. 2(c)], the elastic centroid the plate the and directions are 
not located the same plane parallel the face plates. Asa result, the formulation 
for such plates are derived with respect reference plane, (Fig. 1), which 


EM2 
___ 
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does not contain the elastic centroids all directions. The coupling stiffnesses, 
[C] (1), accounts for the difference locations the elastic centroids 
the plate and the difference plane chosen. Coupling stiffness, [B,] (1), also 
exists when the reference plane face plate, (Fig. 1), does not contain 
the elastic centroids the face plate the and the directions [see Fig. 
L1, L2, and L3, and Fig. For unbalanced cross-ply laminated faces 
[Fig. the resulting coupling stiffnesses not alter the orthotropic nature 
the sandwich plate. However, for unbalanced angle-ply laminated faces [Fig. 
the resulting coupling stiffnesses are anisotropic, which has made solutions 
for such plates inherently abstruse. Note that neglecting the anisotropic coupling 
stiffnesses, and BY? (1), can result appreciable errors the deflection 
and frequency predictions, the magnitude which depends several factors 
related the elastic properties the face plates and the core and the geometry 
the sandwich plate. 

order establish the magnitude the coupling effect the deflection 
and natural frequencies unbalanced angle-ply sandwich plates, the center 
deflection, the sandwich plate subjected uniform transverse load 
in. (6.9 and the natural frequency, w,,, the plate are 
compared with the results for the center deflection, w,, and the natural frequency, 


respectively, equivalent sandwich plate with the anisotropic coupling 


FIG. 3.—Deflection Ratio, and Frequency Ratio, Function Thickness Ratio, 
and Ply Orientation, for GR-WC Plate 2.5) 
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FIG. 4.—Effect Face Stiffness Deflection Ratio, and Frequency Ratio, 


FIG. Ratio, Function Aspect Ratio, Thickness Ratio, 
and Ply Orientation, for GR-WC Plate 
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this study, the face plates are constructed either graphite-epoxy (GR) 
glass-epoxy (GL) angle-ply composites, having the properties listed Table 
The cores considered include weak core (WC) and relatively stiff core 
(SC), with the properties given Table all the cases considered, the 
core thickness taken in. (25.4 mm), while the thickness the face 
plates vary from 0.01 in.-0.25 in. (0.254 mm-6.35 mm). the case thin 
angle-ply plates (5), the influence bending-membrane coupling diminishes 
rapidly with increase the number plies. Therefore, order generate 


extreme coupling, each face considered consist only two equal thickness 
plies (see Fig. 2). 


FIG. 6.—Frequency Ratio, Function Aspect Ratio, Thickness Ratio, 
and Ply Orientation, for GR-WC Plate 


The effect the ratio the thickness the core that the faces, 
demonstrated Fig. for GR-WC plate with aspect ratio 
2.5, for various ply orientation, can realized that the thickness 
ratio, has major influence both the deflection and frequency ratios. 
Inspection Fig. reveals that the effect the anisotropic coupling stiffness 
both the deflection and the natural frequency the plate practically vanish 
for very thin faces 20). the other hand, for relatively thick faces 
10), the coupling effect appreciable. 

pointed out that coupling always reduces the stiffness the angle-ply 
sandwich plate considered herein. That is, neglecting the anisotropic coupling 
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stiffnesses can generate errors the deflection and frequency results which 
can reach high 150% for the deflection and 50% for the frequency predictions 
(Figs. 3-7). general, the effect coupling the deflection response 
always higher than its effect the frequency response angle-ply sandwich 
plates. 

Note that, since the magnitude the anisotropic coupling stiffnesses increase 
with increase the face elastic moduli well the ply orientation, 
the effect coupling the deflection and the frequency, always increases 
with increase the elastic moduli (Fig. and the angle 45) 
(Figs. and 7). 


FIG. 7.—Effect Core Stiffness Deflection Ratio, and Frequency Ratio, 


The effect the aspect ratio the plate, L/B, the deflection 
ratio, and the frequency ratio, are presented Figs. and respectively. 
can realized that the effect anisotropic coupling stiffnesses both 
the deflection and the frequency responses increase with increase the 
aspect ratio. 

The effects the core shear stiffness the deflection and the frequency 
are presented Fig. from which can noted that the coupling 
effect increased with decrease the core shear stiffness. 

general, can stated that the effect the anisotropic coupling stiffnesses 
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the deflection and the frequency responses depends mainly upon the relative 
thickness the core and the faces (7), and that the effect magnified with 
increase the aspect ratio, the stiffness the face plates and ply orientation, 
and/or decrease the core shear stiffness. Numerous numerical results 
this study verify that the effect shown Figs. 3-7 typical. 

order quantify the previous analysis, results for the natural frequencies 
sandwich plates with unbalanced cross-ply and angle-ply faces are presented 


TABLE 3.—Natural Frequencies for Cross-Ply Sandwich Plates in. (127.0 


Model Number Lay-Up Design 


TABLE 4.—Natural Frequencies for Angle-Ply Sandwich Plates in. (127.0 
cm), 1/8 in. (0.318 cm), 1.0 in. (2.54 cm)] 


Model Number 


— 


Tables and respectively. The core considered orthotropic with shear 
cross-ply cases, solutions are generated for the ply lay-ups L2, and 
shown Fig. The results presented Table are excellent agreement 
with those obtained using the formulations presented Refs. and For 
the angle-ply cases, the results obtained were not verified, since the writer’s 
knowledge, rigorous solutions are unavailable for such plates. expected, 


(3) (4) (5) 
152.6 146.0 159.3 
422.7 401.1 444.3 
670.8 744.1 
197.6 196.4 199.2 
442.4 423.7 461.5 
720.0 685.3 754.9 
273.7 279.5 268.7 
482.0 468.8 496.2 
745.1 714.2 776.7 
Angle 
(5) 
90.60 122.6 152.9 
199.9 283.3 382.1 
371.0 519.5 658.7 
179.4 209.2 227.3 
275.0 422.9 
423.2 542.1 678.4 
276.1 300.4 308.6 
365.2 428.7 484.4 
499.7 591.1 711.6 
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however, that the degree accuracy the results the angle-ply cases 
the same those for the cross-ply cases. 


The modified stiffness formulation extended the dynamic analysis 
sandwich plates with anisotropic unbalanced laminated faces, such those 
with unbalanced angle-ply faces, which were not amenable rigorous solutions. 
The formulation independent the boundary conditions and can used 
conjunction with variety analysis techniques, such the finite element 
method. Using double Fourier series solution, the method has been applied 
the analysis sandwich plates with unbalanced cross-ply and unbalanced 
angle-ply faces. The results for the cross-ply cases are complete agreement 
with those obtained using other analysis techniques. expected, therefore, 
that the solutions for the angle-ply cases have the same degree accuracy 
those with cross-ply faces. 

was pointed out that the difficulty associated with the analysis sandwich 
plates with unbalanced angle-ply faces due the existence anisotropic 
coupling stiffness the formulation. The study the effect anisotropic 
coupling stiffnesses the deflection and natural frequency sandwich plates 
with unbalanced angle-ply faces reveals that for certain types sandwich plates, 
the effect coupling significant. For such plates, neglecting coupling can 
result appreciable errors which can reach high 150% for the deflection 
results, and 50% the frequency prediction. can concluded, general, 
that the overall influence the anisotropic coupling stiffnesses the behavior 
unbalanced angle-ply sandwich plates governed mainly the thickness 
ratio was found that for relatively thin faces 20), the coupling 
phenomenon practically vanishes; however for relatively thick faces 10), 
the coupling effect appreciable and was found that the effect coupling 
magnified with increase the aspect ratio the plate, the stiffness 
the faces, and the angle orientation the plies, and/or decrease 
core shear stiffnesses. Finally, since the effects factors which influence 
coupling are interdependent, the anisotropic coupling stiffnesses should always 
considered the analysis sandwich plates with unbalanced angle-ply 
faces. 
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The following symbols are used this paper: 


x? 


M”, MY,M2 


™ sy 


Pe» Py 


Superscripts and Subscripts 


surface area plate; 

matrices containing sum stiffnesses; 
membrane, coupling, and bending stiff- 
nesses; 

width plate; 

modified stiffness matrices; 

shear moduli orthotropic core; 

inertia constants associated with rotary iner- 
tia core and faces, respectively; 

length plate; 

lay-up design; 

vector containing sum moment resultants 
two-face plate; 

moment resultants face plate, 

inertia constants associated with translatory 
inertia core and faces, respectively; 
transverse shear forces core; 

total kinetic energy; 

kinetic energies core and faces, respec- 
tively; 

thickness core and face plates, respec- 
tively; 

strain energy two plates; 

total transverse deflection; 

partial transverse deflection and direc- 
tions due bending, respectively; 

partial transverse deflection and direc- 
tions due shear, respectively; 
undetermined displacement coefficients; 
rectangular coordinate system; 

angle between principal axes lamina and 
reference axes; 

vector partial curvatures; 

densities core and faces, respectively; and 
frequencies. 


core sandwich plate; and 
face sandwich plate. 
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TECHNICAL NOTES 


BUCKLING RINGS SUBJECTED 
INTERACTING 


INTRODUCTION 


This technical note deals with the interaction effects multiple external 
and internal pressure loads, e.g., normal gas pressure, centrally directed pressure, 
etc., acting circular elastic ring buckling its own plane. Several simple 
formulas for the critical load are derived and number unusual phenomena 
are discovered. Moreover, misconception concerning the buckling ring 
subjected the constant-directional pressure load clarified. 

Heretofore, the extensive technical literature the stability rings has 
dealt only with single distributed pressure load multiple point loads. 


possible imagine several different ways applying fairly uniform 
external internal pressure load ring. The most common way—at least 
the published literature—is subject the ring gas hydrostatic pressure 
that remains normal the ring buckles. Pressure loads with non-normal 
directions during buckling can produced attaching strings the ring 
closely spaced equidistant points; placing the ring horizontally smooth 
table; drilling concentric circle holes, corresponding number the strings, 
through the table; passing the strings through these holes; and then attaching 
weights the free ends the strings, thus obtaining essentially conservative 
load system. drilling more than one circle holes, and passing strings 
through each them, multiple interacting pressure load systems can obtained. 

Herein, order save space, all steps the derivation the governing 
equation will omitted, these are elementary, and can found many 
publications, e.g., Refs. The derived governing equation 


i=! 


i 
which circumferential (tangential) displacement point the centroidal 
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line; radius the circular centroidal line the ring; 
being the position angle measured from any reference radius, 
elasticity, centroidal moment inertia the cross-sectional area, 
external gas (normal) pressure load per unit length, and inward (initially 
radial and normal) pressure load exerted the strings passing through given 
circle holes; distance measured inward along diameter from the ring’s 
centroidal line point the given circle holes; and indicates summation 
all pressure loads produced the strings passing through the holes 
the several circles holes (each will characterized different direction 
after buckling). 

Since the solution Eq. for closed ring must periodic, posses 
certain symmetry, and devoid rigid-body displacements 


Substitution Eq. Eq. results the algebraic equation 


a 
i 


With the notation 


N 


Eq. can rewritten the form 


which can easily solved for the dimensionless critical pressure load 
terms the parameters and R,. 

when the ring subjected external pressure always directed 
The three different critical values calculated agree with those frequently 
quoted the technical literature, e.g., Refs. 5-7. 


Consider circular ring subjected two interacting pressures and p,, 
Then, from Eq. for 


a 
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Equation indicates that the ring under this combination loads can 
also observed that for the dimensionless critical net pressure 
load positive (compression) and decreases with increasing and k,. (The 
stability boundaries consist straight line segments and not violate any 
the known theorems their properties.) 


Consider circular ring subjected two interacting pressure loads and 


constant-directional pressure and centrally directed (polar) pressure. 
Then, from Eq. for 


Equation indicates that the ring subjected this combination loads can 


Consider circular ring subjected only one pressure load, say 
exerted the strings passing through the holes the circle holes characterized 


and yields 108/11 for and exactly the same value 108/11 
for This result implies that for 3a/16, such string-loaded ring 
will buckle higher mode than the first. (It should noted this point 
that, for very small values the ratio a/h,, the ring’s buckling behavior will 
similar that shallow arches for which the assumptions made this 
communication are invalid.) 


Remarks 


The foregoing three problems obviously not exhaust all the possibilities 
combining different pressure loads ring; they merely serve illustrative 
examples. 

Now, the ring that subjected the constant-directional pressure not 
free, but hinged two diametrally opposite points, its critical buckling pressure 
Ref. However, even this critical value, obtained via the general solution 
Eq. and the appropriate boundary conditions, not the lowest. careful 
analysis hingeless circular arches means Eq. indicates that when 
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the ring fixed one point against both rotation and translation, 0.7014, 
the case constant-directional pressure—a truly surprising result. 
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SOLUTION 
INVERSE PROBLEMS 


INTRODUCTION 


Many problems engineering and science are converted system linear 
simultaneous equations, i.e.: 


[A] {x} 
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Such set equations may arise approximation continuous relationship, 
expansion the continuous functions y(t) and terms appropriate 
sets orthogonal functions, which case and represent expansion 
coefficients. 


the functions are linear x,, may write the problem the 
form 


the matrix form Eq. the functions are not strictly linear, 
but vary smoothly enough, they may expanded Taylor series about 


some set initial values the say 
o . 


x? 


Defining Ay, and ignoring second and higher order terms 
have 


This the same form Eq. with the substitution Ay, for y,; Ax, for 
x,; and 


For simplicity, shall proceed using the notation Eq. with the understanding 
that the aforementioned substitution can made any stage the calculations 
for system which results from the perturbation quasilinear problem (3). 

Eq. well-posed problem, the solution this system simultaneous 
equations obtained any the wellknown existing techniques (1). However, 
many circumstances Eq. cannot solved directly. This occurs when matrix 
ill conditioned rectangular matrix. There are methods available 
handle this situation, such Orthogonal Decomposition Methods (4), and 
Natural Inverse Matrix (5); but any event, the solution obtained for 
any direct method often not acceptable. This due the fact 
that the error introduced the conversion from continuous system discrete 
system has not been considered the solution. 

The objective this paper present method solution which will 
account for the discretization errors the solution. The technique will 
compared the present method solution which has been used many 
cases for acceptable solution Eq. 

Consider the following optimization problem: 


[c] {x} 


0A, 
Ax, 
Ox, x? 
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and {e} {e} 


which [c] banded positive definite square matrix, varies from problem 
problem, and the limit might identity matrix; and {e} vector 
upper bound acceptable error. 

Equation may obtained from continuous minimization problem after 
discretization. The Lagrangian function Eq. 


and with the gradient this Lagrangian function set zero 
VL(x,y) 


which y~' Lagrangian multiplier. the Lagrangian Multiplier y~' known, 


substituting into Eq. and solving for {x}, one will obtain the optimal 
solution the optimization problem (Eq. 4). However, the Lagrangian multiplier 
not known, and the known quantities are the upper bound acceptable 

Therefore, there need for method solution the primal optimization 
problems (Eq. 4). The following algorithm technique for finding optimal 
solution Eq. This algorithm has the advantage solving convex nonlinear 
optimization problem series linear programming optimizations. 


can shown that the optimization problem (Eq. convex minimization 
problem, (2,6). Any convex minimization problem has unique solution. 
solution algorithm solve convex minimization problems follows. Consider 
the general convex minimization problem: 


Minimize 
Subjectto 
which f(x) convex function; g(x) vectors concave functions; 
and set all feasible solutions: 
Construct linear programming problem: 
Minimize 


Subjectto 


Solve this linear programming problem for primal solution and dual solution 
Construct minimization problem: 


Minimize f(x) 
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Subjectto 
Let say the optimal solution Eq. 10. 


problem (Eq. 9). 


Minimize 


Let say that are primal and dual optimal solutions Eq. 


Construct minimization problem 
Minimize f(x) 
Subjectto 


with new Lagrangian multiplier and solve the problem. Let say 
the answer this minimization problem. 

programming problem (Eq. 11) and continue. For proof and convergence 
algorithm see Refs. and This algorithm will give the optimal solutions 
the primal problem (Eq. and the dual system (Eq. 10) simultaneously. 


There are several important technical cases where the objective solve 


and least square error procedure was used. 


order obtain acceptable smooth solution, positive definite banded 
symmetric matrix added the matrix such that 


which [c] banded positive definite symmetric matrix; and small 
positive number. applying this technique the intention include some 
error the solution, but the result not always satisfactory. this method 
not possible control the contribution the error the problem 
direct way, trial and error method should used obtain the desired solution. 
course the solution Eq. depends parameter and varying 
one would have set solutions, which may may not feasible solutions 
and none them are the optimal solution Eq. 


The two-story frame structure Fig. with element properties 
A(t) 100 sin 
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A(t) 


FIG. 1.—Two-Story Frame Structure Subject Harmonic Ground Excitation 


1.00 


Normalized Unit Impulse Function 


FIG. 2.—Unit Impulse Response Function: (a) Exact Solution; Direct Solution; 
and (c) Optimal Solution Integral Eq. 
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The absolute acceleration y(t) the roof level related the ground 
acceleration 


0.05 and integral Eq. discretized: 


j= 


which depends the numerical integration scheme used discretization. 
this example, 40-points time increment was considered. Equation resulted 
simultaneous equations shown Fig. 2(b), where the exact answer 
shown Figure 2(a). The optimization technique applied this problem and 
further details each step are given Ref. The optimal solution shown 
the Fig. 2(c). 


method solution and algorithm were presented this paper for the 
solution many converted problems engineering and science system 
linear simultaneous equations. The method based minimizing objective 
function subjected the upper bound the error the constrained equations. 
The optimization problem one convex minimization, and unique solution 
the problem exists. The method applied two-story frame structure 
subjected harmonic ground excitation and the optimal solution good 
agreement the exact solution. 
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the Publications ASCE for more detailed information preparation and submission 
manuscripts. 


DISCUSSION 


The writer appreciates the comments Stathopoulos and will endeavor 
clarify our presentation the results the study, hopefully his satisfaction. 
and other readers should refer once the section this 
Closure for correction several numerical errors which unfortunately appeared 
the original paper for which the writer sincerely apologetic. 

The flows discussed are parallel the face the building, indicated 
the Introduction, and are treated terms two components: one parallel 
the vertical mullions; and the other horizontal. The component normal 
the building surface and away from mullion influence zero. Thus, the 
flow 90° the mullions, purely horizontal flow parallel the 
building face while, the flow 45° the mullions, its vertical and horizontal 
components are equal magnitude. The appropriate model orientation then 
with plate parallel the tunnel center line and mullions either transverse 
45° the tunnel center line. actuality the plates used were plywood 
above the tunnel floor and spanning its 3-ft (91.4-cm) width with streamwise 
slope adjusted for constant static pressure along the tunnel roof. Upon one 
plate were mounted six 1/4-scale I-beams in. (5.08 cm) high and set 90° 
the tunnel center line. Upon the second plate were mounted six 1/4-scale 
mullions but 45° the tunnel center line. Upon third plate were mounted 
1/10-scale mullions 45° the center line. 

For the 1/4-scale, 45° model the ratio mullion spacing, tunnel 
width, was 0.417. This value L/W was shown too large for free 
lateral adjustment the flow the disturbing influence the mullions 
the observation that the pressures measured the plate were nearly the same 
those which had been measured the 90° model. 

For the 1/10-scale, 45° model the value L/W was 0.167, and the 
differential pressures measured were less than half those measured the 
1/4-scale, 90° model. This was attributed the freedom the streamlines 
near the plate deflect the direction the mullions and exert pressure 
the mullions only proportion the square the normal component 
velocity. When pressure coefficients were calculated this way, agreement 
was within 0.03 all points excepting that 1.0. not reasonable 
speak the discrepancy relative local pressure coefficients which range 

1979, Bernard Leadon and Michael Kownacki (Proc. Paper 14732). 


Engrg. Sci., Univ. Florida, 231 Aerospace Engrg. Bldg., Gainesville, Fla. 
32611. 


Engr., Goodyear Atomic, Piketon, Ohio. 
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from negative positive unless possibly comparison with the whole range 
which 0.26; thus the discrepancy only 11% except, again 1.0, 
where 35%. Put another way, the mean slope the versus distribution 
for the 1/4-scale, 90° model was approx 0.27 over from 0-0.917. 
For the 45° model was approx 0.24 when the normal component 
velocity was used compute the dynamic pressure. This slope would have 
been approx 0.10 the full free stream dynamic pressure were used. these 
terms then stated that the distributions for the 1/10-scale, 45° 
model fairly well’’ with the distribution for the 
90° model, and ‘‘at 0.167 the independence principle appears 
extended study, course, one might test additional scales model 
ascertain with greater precision the effect upon measured pressures. 

From the results just discussed, however, the ratio L/W clearly shown 
important parameter wind tunnel testing two-dimensional planar 
geometries, much for three-dimensional bluff bodies the model tunnel 
cross-sectional area ratio important parameter. either case the parameter 
only indicates maximum value which must must not exceeded the test 
similarity parameter which must held constant between model and prototype. 
the other hand taken for granted that the geometries model and 
prototype, including the ratio mullion spacing mullion height, must 
identical for geometrical similarity. 

the numbers mullions the two models, filled the available 
model length with mullions which resulted twice many for the smaller 
scale model. There reason hold this number constant since after four 
mullions the turbulent boundary layer and separated flow zones between the 
next two will very little affected adding more mullions forward. 
testing smaller model proportionately smaller tunnel, one would certainly 
expect duplicate results already obtained both tests were done with 
comparable accuracy. 

Data extreme pressure coefficients were not recorded, but agreed 
that such measurements could interesting. 

With regard the effects turbulence the free stream were convinced 
that four more upstream mullions would generate enough turbulence the 
boundary layer render mean surface pressures independent free stream 
conditions, but were perhaps negligent not explaining this point. has 
long been known, e.g., that heat transfer turbulent boundary layer 
plate exposed nonaccelerating freestream quite insensitive freestream 
turbulence level (3). have recently had occasion verify this point and 
have reported elsewhere (4) that the change from 0.05% 4.5% turbulence 
level makes difference whatsoever the pressure coefficients measured 
the 1/4-scale, 90° model. 

are grateful for the opportunity present this closure. also note 
that the National Science Foundation has supported this work under Grant 
No. PFR-7726391. 


Kestin, J., Maeder, F., and Wang, E., Turbulence the Transfer 
Heat from Plates with and without Pressure Journal 
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Heat and Mass Transfer, Vol. No. Sept., 1961, pp. 133-154. 
Leadon, and Jenkins, A., Curtain Wall, Full Scale and 


Model Results,’’ presented the October 27-30, 1980, ASCE Fall Convention and 
Exhibit, held Hollywood, Fla. 
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Errata.—The following corrections should made the original paper: 


Page 515, paragraph line Should read ‘‘may modeled wind 

Page 517, paragraph line Should read ‘‘The tap locations were x/L 
0.038, 0.25, 0.417, 0.583, 0.75, 0.917, and instead tap 
locations were 0.083, 0.25, 0.417, 0.50, 0.583, 0.75, 0.917, and 

Page 517, Fig. 2(b): delete points plotted 0.083 and 0.417 for 
90° 1/4-scale model. These points are correctly plotted Fig. 2(a) and Fig. 
3(a). 

Page 519, paragraph line Should read calculated was 

Page 519, paragraph line Should read one case was 90° 
while the other was instead one case was zero 
while the other was 

Page 520, line Should read for the 1/4-scale instead 
for the 1/10-scale 


The writers would like thank Bazant for his interest in, and discussion 
of, their paper. While they agree that orthotropic models may eventually 
shown inadequate represent the complexities concrete behavior, 
they believe that the arguments advanced Bazant basis for rejecting 
all orthotropic and hypoelastic models are fallacious. 

Bazant’s arguments appear arise out lack appreciation number 
fundamental characteristics the nonlinear hypoelastic orthotropic models 
that are currently being used simulate concrete behavior. Some these 
characteristics may itemized as: 


For any orthotropic relationship normal and shearing stresses and strains 
are coupled except the principal directions orthotropy. 


1979, Alaa Elwi and David Murray (Proc. Paper 14746). 


Engr., Civ. and Structural Dept., associated with Pullman-Kellog Ltd., 
Edmonton, Alberta, Canada, T6G 2G7. 


Civ. Engrg., Univ. Alberta, Edmonton, Alberta, Canada, T6G 2G7. 
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Orthotropic behavior assumed only incremental sense which does 
not imply orthotropic relationship between total stresses and total strains. 

The principal directions orthotropy are assumed rotate loading 
progresses. The direction incremental orthotropy independent assumption 
the model. 


The arguments and illustrative example presented Bazant demonstrate 
lack invariance are fallacious because they not recognize these characteris- 
tics. Given any assumed directions incremental orthotropy, increments 
stresses and strains and the incremental constitutive coefficients are tensor 
quantities and there lack invariance. one can categorically reject 
all incremental hypoelastic relationships how can one then accept incremental 
plasticity relationships which are special cases (irreversible) hypoelastic 
relationships (20,21). 

The writers have arbitrarily chosen the directions incremental orthotropy 
the principal directions total stress. This appears work reasonably 
well for the test data available them. Whether not this assumption will 
prove acceptable simulate data from tests which the principal axes rotate 
can only verified after reliable data this type are experimentally produced. 
the meantime the writers maintain that choosing the principal stress directions 
the directions orthotropy imposes condition the material law 
which, rather than limiting the capability the model, stabilizes it. There 
implied coincidence between the principal directions total stresses and 
total strains nor between the principal directions incremental stresses and 
incremental strains. general finite element formulation the increments 
equivalent uniaxial strains are accumulated but the increments stress are 
not. Rather stress point which satisfies the stress-equivalent uniaxial strain 
curve found, and equilibrium established iterative manner. The 
instantaneous incremental constitutive matrix depends upon the entire history 
the material point, since function the equivalent uniaxial strains. 
Therefore induced damage affects the current well future values strains 
and also the orientation principal stresses. 

his fourth paragraph Bazant makes number statements regarding 
nonvanishing shear stresses and isotropy. The writers cannot accept these 
statements. They believe has misinterpreted this model, well other 
hypoelastic models (2,4,5). 

The writers also believe that dilatancy has little direct structural significance. 
When dilatancy suppressed confinement, the resulting behavior can 
observed directly biaxial and triaxial tests and this behavior represented 
the stress-equivalent uniaxial strain curve, the ultimate surfaces, and 
the representation Poisson’s ratio. 

The writers agree with Bazant that more basic experimentation including shear 
components necessary. However the writers believe that results from biaxial 
and triaxial tests form the only reliable data upon which base general constitutive 
models the present time. 

The writers also wish point out that endochronic and plasticity theories, 
regardless their success representing concrete behavior, not preclude 
other lines investigation. 


Finally the writers wish admit inaccuracy referring their type 
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ultimate surfaces surfaces. These surfaces would more 
properly referred surfaces (17). 


19. Elwi, E., and Murray, W., Analysis Axisymmetric Reinforced 
Concrete Structural Engineering Report No. 87, University Alberta, 
May, 1980. 

20. Eringen, C., Nonlinear Theory Continuous Media, McGraw-Hill Book Co., 
Inc., New York, N.Y., 1962. 

21. Green, E., and Plasticity Journal Rational Mechanics and 
Analysis, Vol. 1956, pp. 725-734. 


Errata.—The following corrections should made the original paper: 


Page 624, Eq. Should read 


2 


instead 


symmetrical 


Page 637, paragraph line Should read ‘‘and Willam, al. instead 
‘‘and William, al. 

Page 638, Ref. Should read H., Krempl, E., and Willam, 
instead H., Krempl, E., and William, 

Page 639, Ref. 17: Should read ‘‘Willam, J., and Warnke, instead 


symmetrical 
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The importance expedient procedures determine the probability distribu- 
tion the response lightly damped systems subjected nonstationary random 
excitation unquestionable. this context, Spanos and Lutes present their 
paper very interesting results, which may become valuable applications 
random vibration methods many engineering problems. 

The writer believes, however, that some clarifications are necessary. Moreover, 
comments Spanos and Lutes some suggestions presented herein may 
enhance the value their results. 

first place, appears that Eqs. and hold only approximately, provided 
that both a(t) and are slowly varying functions time. Since such assumption 
introduced later the paper, the question remains whether that indeed 
the case. 

Also note that S(t,w,) indicated Eq. 21, the process 
white noise process, and not more general type stationary process. 
such case, the variance the response approximately given 23. 
The writer showed Ref. that the excitation amplitude modulated 
white noise with spectral density function S,, then the variance the response 
x(t) given by: 


which f(t) the modulating function; and When f(s) 
unit step function the origin, Eq. results: 


sin 


which, upon assuming that reduces Eq. 23. 
Formally, for slowly varying modulating function, Eq. can directly 
obtained from Eq. dropping the term cos 2w,(t s)]. This suggests 
1980, Polichronis-Thomas Spanos and Loren Lutes (Proc. Paper 15313). 


*Prof. Civ. Engrg., Curso Engenharia Civ., Universidade 
Federal Rio Grande Sul, Porto Alegre, Brasil. 
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that the expression: 


may more accuracte than Eq. 29. fact, Eq. was valid, the assumptions 
that and that slowly varying function time might relaxed. 
must kept mind, however, that the results apply modulated white 
noise, and therefore constitute only approximations for other nonstationary 
excitations. 

Finally, should pointed out that there error Eq. 29, which 
the argument the factor exp negative. 


25. Riera, D., Método das Estado Estudo Sistemas 
sob Aleatérias Paper No. A-22, Anais 
Congresso Brasileiro Eng. Brasil, 1977, pp. 255-268. 


One the major problems concrete construction determine when 
the in-place concrete has reached sufficient strength that the formwork and 
shoring can safely removed. the past, contractors have used 
maturity aid them. When the weather warm, the formwork and 
shorings are removed early; when the weather cool, they are removed later 
than normal. 

Although this field method similar the maturity concept, some differences 
are obvious. The method just described uses the ambient temperature surrounding 
the forms. Casual acknowledgment the temperatures cannot provide accurate 
data. The definitions warm, cool, early, and late are also left the individual 
contractor. Thus, the contractors remove the formwork and shoring when they 
feel that they can safely so. The removals are not based minimum 
maturity factor but rule thumb past experiences, both. Normally 
this method adequate, however, several recent failures have been attributed 
early removal formwork and shoring when the in-place concrete was 
insufficient strength. Applications exact maturity method appears 
needed within the construction industry. 


1980, Tarun Naik (Proc. Paper 


Assoc. Prof. Construction, Dept. Engrg. and Tech., Rm. 0332, Science Building, 
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The maturity method presented the author can provide accurate 
estimate the in-place concrete strength provided used correctly. The 
user this method must aware the limitations associated with the method. 
Unlike other nondestructive testing methods mentioned the author, the maturity 
method does not measure any the physical properties the concrete (e.g., 
the Swiss hammer measures surface hardness, the Windsor probe measures 
penetration resistance, pull-out test measures pull-out stress, pulse velocity 
measures the velocity transmitted low frequency ultrasonic sound wave). 
The maturity method presents time-temperature history the in-place concrete. 
This time-temperature history presented maturity number can used 
estimate strength provided expected concrete strength the exact mix 
design known. This shown the author Fig. Each mix design provides 
different maturity-compressive strength curve. When the design strengths 
the mix designs are unknown, the maturity concept would inaccurate. 
concrete construction, misbatched concrete (such concrete with lower cement 
factor than specified the design mix), misplaced concrete (placement 
concrete with design strength 3,000 psi where concrete with design 
strength 5,000 psi was specified) mishandling the concrete (poor vibration, 
segregation, adding excess water, etc.) are not uncommon. these particular 
instances mentioned, the maturity method would indicate higher strength than 
the in-place concrete had achieved. Thus, basing decisions for removal 
formwork shoring only maturity method could result failure. 

safely use the maturity method, some method documenting the design 
strength the mix design the concrete must incorporated into the test 


plan. Use other nondestructive tests along with the maturity method would 
appear hold promise providing accurate estimate in-place concrete 
strength designated time. 

The writer shares the author’s opinion that the maturity method has great 
potential. also the writer’s opinion that the user must realize the limitations 
this method and use the maturity method only when the mix design the 
concrete known can determined other test methods. 


Discussion Owen 


The author has done what too rare today: reported extraordinarily 
precise study methods determining in-place concrete strengths. The issue 
critical for construction safety, the strength little rapid highrise concrete 
officially recognized until 28-day cylinders are broken and several floors 
use. 

The earliest age 4-in. 8-in. cores, drilled days and dried for 
before testing, are said some higher strength (based maturity 
and other test relationships) than the actual core indications. Bloem (9) 
reported similar finding for cast-in-place [American Society Testing and 
Materials (ASTM) Method 873] pushout cylinders. the latter 873 pushout 
molds had been used, compressive tests might have been conducted 


Principal, Owen Richards Testing Equipment, R.D. Box 156, The Plains, Va. 22171. 
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earlier ages—when constructors most need significant strength test data 
aid judgements for curing requirements. 

further advantage the 873 pushouts (aside from cost, small fraction 
that for drilled cores) that they are already the as-cured condition, 
thus needing drying out, for wet-drilled core. Early age cores are especially 
subject aggregate-cement paste bond damage and low strength indications, 
before drilling. 

especial problem with evident difference logical approach 
inherent Naik’s Table and what this writer considers significant 
for evaluating probable in-place concrete performance. The writer consider the 
artificial concrete purely contractual interest between concrete seller and 
buyer: enables one judge whether concrete ingredients, mixing and 
handling the point the buyer’s becoming owner are acceptable. The 
potential strength should not, mind, used calculate novel 
potential strength (as Col. Table based correlation with core strength 
what effect different concrete. 

any such correlation made (and may have significance evaluating 
the buyer’s handling, placement and curing efforts), would mind best 
the direction calculating in-place core strengths potential strength. 

Conducting such in-place core potential strength calculation (which, 
noted, the writer depreciates), one gets regression equation, for mixes 
No. and compressive strength cores potential strength) 
0.853 317; the indicated estimated strengths are naturally lower than 
those developed Col. Table Since they are lower, they may qualify 
more conservative, and therefore desirable, judging potential in-place 
strength performance and interim load-bearing capacity the concrete involved. 

The author has presented most interesting and informative presentation 
formal approach in-place strength evaluation. one somewhat familiar 
with the general outlines his active testing program the University 
Wisconsin Milwaukee, this writer looks forward his publication similar 
work more significant and earlier test ages, especially that based his 
pullout strength (ASTM 900) work. 

technique evaluating in-place strength development has been developed 
this writer the course comparing various methods testing concrete 
for strength, with emphasis the pullout. well known, concrete strength 
tends increase strength when tested various points starting the top 
and moving toward the bottom specimen structure. The writer hypothesizes 
this trend systematic, confirmed the technique; and that the same 
technique helpful evaluating strength gain versus time and maturity 
indications. 

The technique arose out pullout tests 2-in. square, 5-in. high mortar 
strength cement. The writer tested each column five equally spaced 
elevations: the top (designated level T), middle (M), bottom (B), and upper 
levels (U) and lower (L) levels, midway between and and elevations, 
respectively. Straight-line regression calculations pullout strength 
and on, confirmed the obvious: increasing strengths the 
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approach the could not reach infinity; the pullout strength, 
and positive. 

The same applies strength time: measured successive and about 
equally spaced time intervals, concrete strength gain increments taper off, and, 
naturally, peak out. 

return Naik’s report, all but one the lab cured cylinder lots’ 
strengths conformed the hypothesis. The sole exception was the weakest 
mix, No. which the age formula was 291, with 
the 655-psi incremental strength increase the period exceeding 
all others, the next highest gain period (usually the highest) being the earliest, 
days. 

contrast, none the core lots really met the requirements the hypothesis 
within the usual range, where varies from 0.8-0.9. The 
cores mix No. gave f.. 0.47L 2940, paralleled the same mix’s 
0.64L 2,270, with relatively low slope and high gain rate. 

This the pattern which early-age maturity-strength investigations may clarify, 
if, hoped, they are successful. very early ages strength development, 
which 873 pushout cylinders and 900 pullouts may prove most helpful, 
the slope the formula for strength after successive time periods will exceed 
1.0, reverting the 0.6-0.9 range when the later ages days, with 
which are traditionally familiar, are reached. 

The author congratulated his most informative work. Further publication 


strength development during the period several hours week eagerly 
anticipated. 


Bloem, L., Strength Journal the American Concrete 
Institute, National Ready Mixed Concrete Association Report, Vol. 65, No. Mar., 
1968, pp. 176-187. 
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